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SYNOPSIS 


EXPERIMENTAL INVESTIGATIONS OF SADDLE TO SCISSION 
DYNAMICS IN LOW ENERGY FISSION 


RAGHAVA VARMA 
Ph, D. 

Department of Physics 
Indian Institute Of Technology, Kanpur 
January, 1989- 


Nuclear fission is a process in which a heavy nucleus divides 
in two or more nuclei either spontaneously or induced by an; 
energetic projectile- Discovered in 1939, the early interest in: 
the process was mainly to put the largest known energy release in; 
a reaction (approximately 200 MeV) to commercial and military use-; 
The interest in the physics of the process arose mainly because,; 
apart from giant resonance, fission was the only known nuclear | 
reaction involving rather drastic rearrangement of nuclear matter-; 


While the 

development 

of 

heavy- 

-ion accelerators in 

the 

last 

two 

decades , 

has brought 

out 

other 

new phenomena such 

as 

the 

deep 

inelastic 

col 1 isions 

involving 

nuclear collective 

mo t i on 

to 

the 


forefront, the importance of fission studies has not diminished- 
This is primarily because the nuclear fission can be considered as 
an elementary collective process without the uncertainities of 
reaction dynamics as in the case of heavy ion reactions- Further, 
in the case of low energy fission, the nuclear matter is much 



cooler thermodynamically and therefore the dynamics is strongly 
influenced by nuclear shell structure effects and the superfluid 
nature of nuclear matter- Therefore, fission studies, often serve 
as a test bed of models developed in heavy ion induced reaction 
studies - 

It is known that for all nuclei having mass number greater 
than about 100, fission has a positive Q value- However, only the 
heaviest actinides are known to have measurable spontaneous 
fission rates- This is due to the existence of an energy barrier 
for fission- It is known that this energy barrier can be 
understood in terms of shape evolution of the nucleus and the 
associated deformation potential energy surface- In this model, 
fission is a continuous evolution of the nuclear shapes starting 
from the nearly spherical comptound nucleus to two well separated 
fragments at infinity- In the deformation space, two distinct 
milestone configurations along the. fission path can be identified- 
The first one is the fission saddle point- This is the point of 
maximum deformation potential energy which the nucleus must p«ass 
to undergo fission- In this sense, it is a point of no return from 
where the separation into two fragments becomes imminent- The 
second milestone configuration is the scission point at which the 
nuclear interaction between the two fragments vanishes. The 
fragment mass and charge are uniquely determined at this point- 
How are these affected by the saddle point shapes is still an open 
question- The separation of the two fragments, then takes place 
purely as result of coulomb repulsion- While to the first order, 
the coulomb interaction energy at scission point gets converted 



there 


is also 


into fragment kinetic energy, as they separate, 
shape relaxation of the nascent fragments- The deformation energy 
of the fragments is then converted into their excitation energy- 
The excited fragments then cool by the emission of neutrons and 
gamma rays- In a typical experiment it is the asymptotic fragments 
at infinity which reach the detector. The exp>er imen tal observables 
in fission studies are the relative probability of fission, the 
fragment mass, charge, energy and angular distribution and the 
nature of the secondary radiations emitted by the fragments- The 
nature of the collective motion between the formation of the 
compound nucleus to the saddle point, from saddle point to the 
scission point and the separation dynamics beyond the scission 
point, have important consequences on the final observables - 
Conversely, one can learn about the nature of collective dynamics 
along the fission path from these measurement s - 

It is now known that the motion upto the saddle point is 
highly dissipative and in the case of low energy fission it takes 
place in time scales of the order of 10 seconds- On the other 
hand, the nature of collective dynamics between the saddle and the 
scission points is not fully understood- During this phase, the 
nuclear shape changes from a deformed moriO““riucleus to a higly 
necked di-nucleus- The nature of the dissipative force also 
expected to change correspondingly- While it is estimated that the 
time scale involved, is of the order of 10 seconds or less, the 
precise nature of the dynamics is not known- 

De tailed studies of the secondary radiations such as the 



emitted close to scission have given very important information on 
the scission dynamics- For example, it is now well known that the 
nascent fragments at scission are highly deformed and scission 
takes place off axially, resulting in appreciable angular momentum 
of the fragments from spontaneous fission of even— even nuclei with 
zero spin- However, the effect of pairing, deformation and the 
shell structure are not fully understood- 

The measurement of total kinetic energy of the fission 
fragments as a function of the excitation energy of the target, 
yields information on the degree and the nature of dissipation 
present in collective nuclear dynamics- If the system is viscous 
then the excitation energy of the compound nucleus will not appear 
in the kinetic energy of the fragments, whereas if it is in a 
superfluid state all the excitation energy will go in as kinetic 
energy of the fragments- The degree of excitation energy appearing 
in the total kinetic energy can give a good measure of the 
viscosity of the system- 

In the p«resent work we have studied s 
<a> the gamma ray multiplicity distribution as a function of 

i 

fragment charge ratio - 

(b> anisotropy of the gamma rays with repect to fragment 
direction of motion as a function of gamma ray energy- 
<c> the correlations of the total kinetic energy of the 
fragments and the excitation energy of the fissioning 
nucleus- 

in an effort to obtain further information on the nature of the 
motion between the saddle point and the scission point, and about 



the scission process itself- 

PROMPT GAMMA RAY MULTIPLICITY DISTRIBUTION AND GAMMA RAY 
ANISOTROPY WITH RESPECT TO FRAGMENT MOTION IN THE SPONTANEOUS 
FISSION OF ^®^Cf 

In the first experiment, we studied the multiplicity 
distribution of the prompt gamma rays emitted from the fission 
fragments in the spontaneous fission of Cf as a function of 
fragment charge ratio- The average number of prompt gamma rays and 
the width of the multiplicity distribution are measured- The 
fission fragments were detected by a solid state surface barnier 

252 

detector kept at a distance of 3 cm from a Cf^ source in a 
vacuum chamber. A high pdrity germanium detector <HPGe>, kept at 
90* to the direction of fragment motion, was used to detect the K 
X— rays emitted by the fission fragments- The K X— ray energy is a 
signature of the charge of the emitting fragment- Two sodium 
iodide detectors were placed in a plane perpendicular to the HPGe 
detector. Coincidence was demanded between the pulses from the 
fission fragment detector. X-ray detector and one or both of the 
gamma ray detectors. The K X— r^y spectrum was recorded in a 
multichannel analyzer first in coincidence with the fission 
fragments, then in coincidence with the fission fragments and 
either of the gamma rays and finally in coincidence with the 
fission fragments and both the gamma rays- After suitable 
unfolding of the spectra, the data was analysed to get the mean 
multiplicity and the width of the gamma 


multiplicity 



di stribut i on as a function of the charge ratio of the heavy and 
the light fragments-. The mean multiplicity indicated some odd- 
even dependence* The width of the , distribution exhibited strong 
signatures of shell structure, deformation and odd— even nature of 
the fragments- It showed odd— even effects for nuclei having closed 
shell structure i-e- for Z=51 to 55- In the deformed region i-e- 
2=56—60 the odd— even effects were absent- Further, the width is 
much lower for nuclei in the region 2=56—59- 

These results can be qual i tat ively understood as follows: 
angular momentum in a deformed nucleus is carried as collective 
rotation, whereas in the spherical nucleus it is carried as 
aligned single— parti cle spins- Consequently, the deformed nucleus 
deexcite by the E2 rotational cascades- These being transitions 
from fixed energy states, do not contribute to the width in the 
multiplicity distribution- On the other hand, in a spherical 
nucleus the deexcitation has to be statistical in nature- The 
width of the distribution arises either due to the width in 
statistical emission or due to the initial spin distribution of 
the fragments- As observed in this work, this gives rise to a 
smaller width for deformed fragments as compared to the spherical 
fragments- 

In the second experiment, we studied the anisotropy of the 
prompt gamma rays with respect to the direction of fragment motion 
as a function of gamma ray energy- For this purpose, the ^®^Cf 
source was mounted in a vacuum chamber and the fission fragments 
were detected by two solid state surface barrier detectors kept at 
a distance of 3 cm and at 90 to each other- The gamma rays were 



detected by a large sodium iodide detector <5 inch X 5 inch) kept 
60 cm away from the source and along the direction of motion of 
one of the fission fragments- Fission fragment energy, gamma ray 
energy and the time interval between the jfission pulse and the 
gamma ray pulse were recorded event by event in list mode on a 
magnetic tape using a computer based data acquisition system- Time 
of flight technique was used to eliminate neutrons since they are 
strongly correlated with the direction of fragment motion and any 
contamination by them will thereby result in positive anisotropy. 

It was found that in the region of gamma ray energy 
300<E^<700 keV, gamma anisotropy is positive with a maximum of 
1-2610-01 at around 520 keV- For higher energies anisotropy 
decreases but remains positive- For energies less than around 300 
keV the anisotropy is slightly negative- Qual i tai vely , the result 
can be understood in terms of the Strutinsky’s prescription that 
the E2 transitions are expected to enhance the anisotropy, whereas 
El statistical gamma rays are expected to decrease the anisotropy- 

A Monte Carlo program which simulates the statistical neutron 
and gamma ray emission from nuclei carrying appreciable angular- 

momentum and energy has been developed to calculate the initial 

* ) 

spin distribution of the fission fragments- The code explicitly 
includes nuclear shell effects on the level densities- It is 
assumed that the nucleus dissipates angular momentum only after 
reaching the yrast line by the emission of rotational cascade- The 


program 

uses 

the known 

expressions 

for 

the 

decay 

widths 

for 

neutrons 

and 

statistical 

gamma ray 

emission 

to calculate 

the 

average 

multiplicity and the width 

of 

the 

multiplicity 



distribution due to statistical gamma rays- We have calculated the 
spin distribution of fragments by subtracting the width due to El 
transitions from the experimentally measured widths- 

MEASUREMENT OF THE VARIATION OF TOTAL KINETIC ENERGY <TKE> OF THE 
FISSION FRAGMENTS AS A FUNCTION OF EXCITATION ENERGY OF THE 
FISSIONING NUCLEUS IN THE REACTION ^®®U<a,o(’f> 

We have studied the kinetic energy distribution of the 

230 

fission fragments in the fission of U induced by inelastic 
scattering of 60 MeV alpha particles from the Variable Energy 
Cyclotron in Calcutta- Correlated fission fragments were detected 
by two surface barrier detectors placed at approximately 4 cm from 
the target in a typical back to back arrangement- Scattered alpha 
particles were detected and identified by a E— AE detector 
telescope kept at an angle of 30“ to the beam axis- Coincidence 
between the AE detector and the fission detector formed the master 
pulse for a five parameter recording, in list mode on a magnetic 
tape, of the pulses from the two fission detectors, the two 
detectors of the telescope and the time difference between the 
fission and the AE detector- The energies of the two correlated 
fission fragments yield the masses of the fragments and the energy 
of the scattered alpha particles gives the excitation energy of 
the fissioning nucleus- Suitable corrections for the energy loss 
suffered by the fragments "due to the finite thickness of the 
target and the backing material, neutron evaporation and the 
recoil of the fissioning nucleus were applied- 



We have measured the total kinetic energy in the excitation 
energy range of 8-5 MeV to 28-5 MeV- In the region 6-5 MeV to 8-5 
MeV, the data is not reliable due to the contamination of the 
elastic peak from the Carbon in the backing- It was found that in 
the excitation energy range of 8-5 to 12-5 MeV, the total kinetic 
energy decreases with the excitation energy and above 22-5 MeV, it 
shows a slow increase with the excitation energy- In between, the 
TKE remains constant with the increase in excitation energy- We 
have also measured the variation of TKE as a function of 
excitation energies for mass separated fragments- We find that the 
slope of the curve TK^ versus E increases as a function of 

X 

excitation energy for asymmetric masses but it does not show any 
significant variation- 

While the experiment allowed us to study the effect of 
excitation energy on tlie total kinetic energy of the mass 
separated fragments, the analysis has brought out a serious 
limitation of not only the present results, but also of the 
earlier results- The correction in the single fragment kinetic 
energy due to neutron evaporation was made on the basis of 
previous experimental results which give neutron yield as a 
function of excitation energy of the fissioning nucleus and the 
fragment mass assuming that all the neutrons are emitted from the 
fully acclerated fragments- This assumption may be incorrect 
because of multichance fission i-e- fission following the emission 
of neutrons from the fissioning nucleus at higher excitation 
energy. The emission of these neutrons does not change the kinetic 
energy of- the fragments as opposed to those neutrons which are 



dimerisiomi has 


emitted from the fragments in motion- Another new 
been added by a series of recent measurements of neutrons emitted 
in fission following fusion in heavy ion reactions- It has been 
established by the presence of enhanced number of pre-scission 
neutrons that a substantial excitation energy of the compound 
nucleus is dissipated by the emission of these neutrons before the 
scission takes place- Correction for the neuti-'ons evaporated from 
compound nucleus was made on the basis of standard statistical 
model with the known expression for the escape widths for the 
neutrons and fission- It was found that this correction changes 
the slope of the TkC versus E curve appreciably- Since the 

X 

experimental information on the dependence of the number of 

scission neutrons as a functiori of excitation energy of the 

compound nucleus and the mass of the fragments is not available, 

we could not correct our results for their presence- However, we 

* 

have demonstrated that it is of crucial importance . to know this 
quantitv^ in the determination of nature of the behaviour of total 
kinetic energy as a function of excitation energy- 

The thesis contains the details of the measurement, analysis 
and the conclusions drawn from them in four chapters- Chapter I 
gives general information about the subject of fission dynamics 
and its relevance as a unique probe of nuclear macrophysics- 
Chapter II contains the gamma ray measurements- Chapter III 
discusses the total kinetic energy measurements- In Chapter IV, we 
discuss our results and how they have contributed towards a better 
understanding of the process and where they have fallen short - 

The appendices give details of the analysis procedure used 
for the total kinetic energy measuremerits and the listing of the 
statistical model code and the other programs developed by us in 
the analysis of our data- 
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CHAPTER 1 


FISSION - THE UNIQUE PROBE OF NUCLEAR MACROPHYSICS 

£1.11 INTR<M>UCTION 

The discovery of nuclear 
division of a heavy nucleus 
liberating a large amount of energy^, opened up a Pandora’s box for 
the nuclear physicists. While the early interest in the process was 
mainly to put the largest known energy release in a reaction, 
nearly 200 MeV/, to military and commercial uses, the interest in 
the physics of the process arose mainly because of its unique 
collective character. In fact, the phenomenon was so different from 
the other nuclea^r reactions known at that time that it led to the 
study of the fission process almost in isolation of rest of the 
nuclear physics. While the development of heavy ion accelerators, 
has diluted this uniqueness to some extent by bringing out several 
other new reaction channels involving nuclear collective motion to 
the forefront, the importance of fission studies has still not 
diminished, primarily because fission is an elementary collective 
process, without the uncertaini ties of the reaction dynamics which 
are there in the case of heavy ion reactions- Further, in low 
energy fission, the nuclear matter is thermodynamically cooler and 


fisson in 1939*, involving the 
into two comparable parts and 
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therefore, the dynamics is strongly influenced by nuclear shell 
structure effects and the superfluid nature of nuclear matter. 
Hence, fission studies often serve as a test bed of models 
developed in heavy ion induced reaction studies. 

On the basis of simple energetics fission is a positive 
Q— value reaction , for nuclei having mass number more than about 
lOO- However, only the heaviest nuclei have measurable spontaneous 
fission rates. This is due to the existence of an energy barrier 
for fission which increases rather rapidly with decreasing charge 
of the fissioning nucleus* The existence of this energy barrier can 
be understood in terms of a shape evolution of the nucleus and the 
associated deformation potential energy surface- In this model, 
developed by Bohr and Wheeler^ on the basis of an idea provided by 
Meitner and Frisch , fission is a continuous evolution of the 
nuclear shapes starting from the nearly spherical compound nucleus 
to two well separated fragments at infinity. FigCl-l] <taken from 
Ref. 6) shows the potential energy as a function of deformation- 
The figure also shows the deformation energy contours as a function 
of two principal deformation coordinates (3^ and (3^- In the 
deformation space, two distinct mil&stone con/igxtrations along the 
fission path can be identified — fission saddls point and the 
scission points Fig. Cl. 1] shows the saddle point as determined by 
the liquid drop model. As we see, the fission saddle point is the 
point of maximum deformation potential energy which the nucleus 
must pass to undergo fission- It is, in this sense, the point of no 
return when separation , into .two fragments becomes imminent. At the 
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Fig. 1.1. Liquid drop model potential energy contour of a fissioning 
nucleus os a function of quadrupole and hexadecapole 
deformation parameters. Potential energy versus defor- 
mation along the fission direction is also shown. 

(Ref. 6). 
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scission point, the nuclear interaction between the nascent 
fragments vanishes thereby freezing the fragment mass and charge - 
The separation of the two fragments beyond the scission point takes 
place purely as a result of the Coulomb repulsion between them- 
yhile to the first order the Coulomb interaction energy at scission 
gets converted into fragment kinetic energies, there is also a 
shape relaxation of the highly deformed nascent fragments- This 
^results in the deformation energy of the fragments at scission 
appearing as their excitation energy. The excited fragments cool by 
the emission of neutrons and gamma rays- It is the cold fragments 
at infinity which reach the detectors in an experiment. The typical 
experimental observables are- the probability of fission, the 
fragment mass, charge, energy and angular distributions, and the 
nature of secondary radiations emitted during the fission process 
or by the fragments after scission. The nature of the collective 
motion between the compound nucleus stage and the saddle point, 
from the saddle point to the scission point and the separation 
dynamics beyond the scission point have important consequences on 
the final observables. Conversely, one can learn about the nature 
of the multidimensional, deformation potential energy surface, and 
the collective dynamics along the fission path from these 
measurements - 

Fig[1.21 (taken from Ref.7) shows the average time scales 
involved during the various stages of the fission process- It also 
shows at what stage the secondary radiations from the fragments 
come when they deexcite. The idea about the time scales® 


involved 



Time scale (non linear) 



Fig, 1.2. Schematic illustration of different stages of the fission 
process and corresponding time scale. 

(Taken from ref. 7) 
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in the motion of the nucleus from compound nucleus to the saddle 
point comes from the lifetime studies- The number of neutrons 
emitted from the fissioning nucleus during its descent from saddle 
to scission along with some theoretical input gives us the time the 
nucleus takes in its descent from saddle to scission- The 
anisotropy in neutron and gamma emission gives us the time scales 
involved for fragments to attain the maximum velocity- 

[l.ai MACROSCOPIC APPROACH TO NUCLEAR COLLECTIVE MOTION IN 
FISSICW ^ 

The macroscopic approach of the liqvLid. drop model <LDM>, with 
a few cosmetic touches here and there, has been the mainstay in the 
description of the rich variety of phenomenon unravelled by the 
fission process with the passage of time- This is primarily because 
a typical actinide nucleus consisting of about 240 nucleons has 
about 720 degrees of freedom- A formal description of all the 
dynamical properties of this system starting from a many body 
Hamiltonian and a solution of the relevant Schrodinger equations of 
motion is however quite impractical since neither the many body 
Hamiltonian is known nor is it easy to solve the equations of 
motion- 

In discussing fission, our interest is not in knowing the 
individual behaviour of the nucleons microscopically but only in 
their overall collective behaviour- The problem is considerably 
simplified if one treats the nucleus on its way to fission as two 
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charged drops- The maximum number of degrees of freedom required to 
describe the collective motion of two rotating and vibrating nuclei 
is 17- A workable model of fission, should therefore select only 17 
collective degrees of freedom out of the total 3A degrees of 
freedom and treat them separately from the remaining 3A-17 int^rruxl 
degrees of freedom- 

In such a treatment, the dynamics is strongly affected by the 
assumptions regarding the coupling between the collective and the 
internal degrees. The coupling could be \t)oah, strong or 
intermediate between the two~ These weak and the strong coupling 
give rise to to the two principal branches of nuclear models, the 
adiabatic**'*” and the . statistical model**'**- If further 

assumptions are introduced in either branch, of which the principal 
one is the smallness of the surface thickness compared to the 
nuclear radius <ie an expansion in powers of A ), we obtain the 

1.3 1 « 4 > 

non— viscous and the viscous versions of the liquid drop model ’ 
This classification of nuclear models brings out the fact that the 
liquid drop model is not to be regarded as opposed to either the 
adiabatic or the statistical model but on the contrary, it is a 
particularly simple version of these models and could be used as a 
starting point for either- It should however be borne in mind, that 
the fundamental question that from first principles whether one 
would expect the weak or the strong coupling approximation to be 
more nearly correct is yet to be satisfactorily answered- On the 
basis of the experimental results we cannot discard either of the 
two lines of thought- Intutively, one expects that at high 



excitation energies the strong coupling statistical model may be 
valid while at low energies, partly because of the superfluid 
nature of nuclear matter, the situation may be different. 

The precise choice of the collective coordinates is crucial in 
determinig the nature of motion- The coupling between the internal 
and the collective degrees of freedom has assumed added 
significance with the discovery of other modes of large scale 
collective motion*^'^^ and the question of degree of viscosity 
prevalent in large scale collective motion is now of paramount 
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importance- There have been some early attempts to explain the 
mass asymmetry in actinide fission on the basis of this coupling 
but the idea did not receive as much attention as it deserved- It 
is only in the last two decades that terms like viscosity, friction 
etc- have come into widespread use- Unfortunately till now our 
knowledge^ of this parameter remains sketchy and far from complete- 
A complete description of the three stages of the fission 
process requires knowledge of three physical quantities s 
(a> the d&formjatiort potonttcLl ener^ landscape i-e- the total 

potential energy of the nuclear system as a function of 

deformation parameters in a multidimensional space since 
several deformation coordinates are necessary to describe 
nuclear shapes along the fission path, especially when the 
fissioning nucleus is strongly deformed- 
(b> the inertial mass parameter of the fissioning system at all 

deformations in order to determine the fission path and the 

fission barrier penetrability. 



9 


<c> the dainpiTxg of the fission mode i-e- its coupling to other 
degrees of freedom especially during its descent from saddle 
to scission* 

The potential energy as a function of the collective 
coordinates can be obtained, for example, from the liquid drop 
idealization of a sum of surface and Coulomb energies. Similarly, 
the inertia (in general a tensor) can be calculated for 

i 

irrotational flow of the nuclear fluid- The potential energy 
surface and the inertial mass term are interdependent*. The damping 
term is needed in the macroscopic treatment of a many body system 
by a few generalized coordinates where the remaining internal 
degrees act as a heat sink- This can be visualized as the flow of 
energy from the collective coordinates to the internal degrees of 
freedom. 

There exist two approximate formulae for the rate of this flow 
based on simplified nuclear models^®’***- They are basically based 
on the fact whether the mean free path is small or large- 
(i>If the mean free paths are small compared to the overall size of 
the system, flow of energy takes place from collective motion into 
internal degrees due to foo-hod:^ collisions between the nucleons in 
the nuclei, resulting in a dissipative flow similar to ordinary 
viscous fluids. 

(ii) If the mean path between nucleons is large, then there is the 
so called one body dissipation which results from the collisions of 
the nucleons with the moving boundary of the nuclear potential 


wel 1 ■ 
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Two elemetary dissipation formulae for the one-body 

dissipation have been derived to describe the rate of energy flow 

from collective to single particle degrees of freedom - iKo wall 

/armxila^^ for a freely communicating shape and th^ windaw formxilc^^ 

for two weakly communi eating shapes in relative motion- The wall 

formula is relevant for convex nuclear shapes without a neck while 

a window formula is for shapes with a strong constriction as in the 

late stages of fission or the initial stages of a nucleus-nucleus 

collision* A generalization of the wall and the window formulae is 

called for when the two pieces are deformable and weakly 

communicating as in the last stages of the fission process- There 

22 

have been some attempts in this direction recently 

For cold nuclei with relatively long mean free paths for 
nucleon— nucleon collisions, the two body viscosity is small and can 
usually be neglected in comparison to the one body dissipation- The 
dynamical evolution of the system can then be obtained from the 
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classical equations of motion : 

d fdL ■! _ dL d$ 

3T[_dq^J dq^ dq^^ 

« = Q/2 

where q^ are the generalized coordinates and the second term on the 
R-H.S- of the equation represents the dissipative term- 

A qualitative discussion of the fission process does not even 
require seventeen degrees of freedom. An intelligent choice of 
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collective coordinates may reduce it further- Swiatecki ^ shape 

parametri 2 ing the nucleus, postulated that a minimum of three 
degrees of freedom are essential to estimate the barrier height and 
also the dynamcal evolution of the process- His three macroscopic 
coordinates ares 

1- the relative separation variable, 

2- the neck thickness variable, 

3- the mass asymmetric constant- 

in this approach, the macroscopic deformation potential energy 
is calculated as a function of shape parameters- It consists of the 
sum of an electrostatic energy, nuclear surface energy and shell 
and pairing correction energies. The potential energy landscape for 
a typical heavy nucleus exhibits two local minimas separated by a 
saddle point pass —one of the local minimas correspond to the fused 
system while the other corresponds to a configuration of separated 
fragments- The term local minima has been used since, in general, 
the potential energy is not stationary with respect , to the mass 
asymmetry. The saddle point is, therefore, only a conditional 
saddle point, with the physical meaning of a mountain pass only if 
the asymmetry is effectively held fixed- For not too large 
asymmetries (asymmetry less than the Businaro-Gal 1 one critical 
mass asymmetry) and no shell effects, the potential energy 
decreases with decreasing mass asymmetry. It therefore drives the 
system towards symmetry when the asymmetry relaxes and the system 
becomes symmetric- The energy becomes stationary with respect to 
the mass asymmetry degree of freedom and the conditional saddle 
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point becomes a truLo •unconcLitional saddle point and the bounded 
inner region of low energy becomes the compound nucleus region- 
Another interesting aspect of the potential landscape is the 
existence of two distinct valleys misaligned with respect to each 
other***— the fission valley of connected shapes and the two 
fragment valley of disconnected shapes- They are separated by an 
energy ridge over a narrow range of relative separation around the 
scission point- The motion in the LDM valley is highly damped and 
the dissipation is dependent on the collective coordinate- The time 
scales involved for motion in the three degrees of freedom are also 
distinctly different — partly because of the differences in the 
relevant inertial masses and partly because of the differences in 
the nature of the dissipative forces- In particular, the motion in 
the neck degree of freedom is the fastest and consequently, the 
neck rupture takes place almost at constant relative separation- 

We will now briefly discuss how the different experimental 
studies of the fission process lead to detailed information on the 
structure of the potential energy surface and the nature of nuclear 
collective dynamics on this surf ace - 

£1.31 THE E>EFORMATION POTENTIAL ENERGY SURFACE 

One of the earliest and the most extensive pieces of 
experimenital information on the fission, process is related to the 
structure of the deformation potential energy surface upto the 
fission saddle point. This is obtained by measurement of fission 
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probabilities, half lives in spontaneous fission and fission cross- 
sections in particle induced fission- Though the height of the 
fission barrier is a small fraction of the total energy released in 
the process, it plays a crucial role in deciding the fission 
probabilities and the angular distribution of the fragments- The 
experimental fission barrier height exhibits a strong dLecrecxse with 
increasing chjcrrge of the fissioning nuclei- This is in qualitative 
agreement with the predictions of the liquid drop model- However, 
the simple liquid drop model fails to explain the near constancy of 
the fission barrier height and the asymmetric mass distribution in 
the fission of actinide nuclei- This alongwith the more recent 
discovery of long lived fission isomers and widely spaced 
sub— barrier fission resonances has led to the development of the so 
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called macroscopic -microscopic approach for the calculation of 

nuclear deformation potential energies- With a deformation 

dependent shell correction energy superimposed on the smooth liquid 

drop model deformation potential energy, it was shown by 
28 

Strutinsky that the fission barriers are doubl e— humped for 
actinide nuclei- In this approach, for all deformations along the 
fission path the potential energy of the nucleus can be expressed 
as a sum of two terms E <s) and E <s>- The E <s> represents the 

M Tn M 

energy obtained by the LDM, whereas E <s> is the 


shell correction 



Proton number 



fig. Cl. 32 Schematic 


Illustration of The Double Humped barrier- 



Neutron number N 


Fig. 1.3(b). Predicted island of superheavy nuclei around the doubly 
closed shell 2*114, N = 184 . 
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which is derived from the value of single particle level density at 
the Fermi surface- Nilsson model with spheroidal shape was used 
to obtain E <s>- FigCl-3a3 shows the liquid drop potential as a 

m 

fuction of deformation now incorporating Strutinsky’s microscopic 
correction to get the double hump^ed barrier- To digress a bit, one 
of the most important predictions of this approach was the 
existence of supei — heavy elements [see Fig-l-3b3- This was one of 
the motivations for the development of heavy ion accelerators - 
Though we have had limited success in the synthesis of transuranium 
elements, the advent of accelerators has provided other 
manifestations of large scale collective motion- 

It has been shown that while the first barrier is mass 

30 

symmetric, the second barrier is mass asymmetric . -Detection and 
measurement of the rotational gamma rays in the second minimmum has 
provided direct experimental support to the picture of the double 
humped fission barrier^- The deformation potential energy surface 
for actinide nuclei can now be calculated to an accuracy sufficient 
to yield fission barrier heights within 0-5 MeV- There have also 
been some recent measurements of the conditional fission barrier 
heights through studies of highly asymmetric fission which are in 
general agreement with the predictions of the liquid drop model - 
Fission barrier heights of medium mass nuclei continue to defy a 
complete understanding with the liquid drop model over predicting 
the barrier heights by 5—10 MeV- Second order corrections to the 
liquid drop model such as finite range effects, compressibility and 
charge redistribution effects seem to be playing an important role 
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in deciding the fission barrier heights 

Further investigations along these lines are clearly required- 
There is no direct experimental information on the structure of the 
deformation potential energy surface beyond the fission saddle 
point- The calculated slope of the potential energy beyond the 
saddle point depends rather sensitively on the magnitude of the 
nuclear curvature tension because of the highly necked in shapes 
encountered at this stage of the process- The value of the nuclear 
curvature tension is however not known unambiguously at present 
because there is a discrepancy between the theoretical estimates of 
this quantity and the value deduced from nuclear mass fits- The 
determination of slope is important as it contains information on 
the nuclear viscosity present during the descent of the nucleus 
from saddle to scission- Considering that the nature of collective 
dynamics beyond the fission saddle point will be dependent on the 
slope of the potential energy, there is an urgent need to resolve 
this discrepancy- However, there are a number of experimental 
measurements which can throw some light on this problem- 

One of these is the study of the light charged particles 
emitted in the so called ternary /is^ion - The angular 
distribution of the LCPs is sharply peaked perpendicular to the 
fission fragment direction- This implies that LCPs are emitted from 
the vanishing neck region between the two nascent fission fragments 
when they are fairly close to each other- Its subsequent motion can 
be understood in terms of the effect of the Coulomb field of the 
two fragments- Based on this a trajectory calculation can be 
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performed which alongwith other things gives us the idea of 
prescission kinetic energy i-e- the energy accumulated by the 
nucleus before it divides in two- The value of the prescission 
kinetic energy is an important parameter in the determination of 
the viscosity- 

Recent measurements on the emission of neutrons in 
coincidence with fission have shown that an appreciable fraction of 
the neutrons are emitted during the saddle to scission descent of 
the nucleus- The fraction of excitation energy going into the 
emission of these scission neutrons as compared to the energy 
appearing in the prescission kinetic energy can yield information 
on the nuclear viscosity during this stage of the fission process- 
However, our present understanding of the reaction dynamics and 
also the experimental input is still not adequate enough to deduce 
a reliable value of the viscosity coefficient- 

A major fraction of the energy released in fission appears as 
the total kinetic energy of the fission fragments- As mentioned 
earlier, this energy first appears as the Coulomb interaction 
energy at scission and then gradually gets converted into fragment 
kinetic energies as they separate to infinity- The variation of 
total kinetic energy as a function of Z is well documented- It has 
been found that to the first order total kinetic energy varies 

“ 2 J-/3 3*4 

linearly with Z /A - Substantial experimental information on 
the variation of total kinetic energy with the excitation energy 
for the actinide nuclei exists in literature®^ -The variation seems 
to depend crucially on the relative height of the two barriers- 
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However, in light of recent results on increased scission neutron 
multiplicity in both heavy ion as well as light ion induced 

reaction , the results on the variation of the total kinetic 
energy with excitation energy of fissioning nuclei have become 
suspect. We will elaborate more on this in Chapter 3 where we 
present our result on the measurement of the variation of total 
kinetic energy with excitation energy for fission fragments from 

238 

U- The nature of this variation determines the magnitude of 
viscosity. The nature of the collective dynamics between the saddle 
and the scission points and temperature dependence of the scission 
act may yield a small dependence of the kinetic energy on the 
excitation energy of the compound nucleus. Experimental information 
on this aspect of the process is quite meagre and additional 
experiments may be most useful- 

£1.41 MASS AND CHARGE RELAXATION IN FISSION 

The asymmetric fragment mass and charge distributions in low 
energy fission of actinide nuclei , have remained a challenge to 
nuclear theorists for several decades- While it has been 
established well that the asymmetry arises mainly due to the shell 
structure effects of the nascent fragments, the exact mechanism of 
mass and charge evolution leading to the observed asymmetry is not 
yet known. While the quasistationary models have had some success, 
it has been established that simple considerations based on the 
static deformation potential energy are not adequate to explain the 
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asymmetry. The collective dynamics in the mass asymmetry degree of 
freedom seems to be playing a crucial role beyond the fission 
saddle point- This includes not only the kinetics governed by the 
inertia and the driving potential but also the nature of the 
dissipative forces in operation between the saddle point and the 
scission point- This has given rise to statistical models based on 
full equilibrium at the scission point, stochastic models of 
nucleon exchange between , the nascent fragments^^ and dynamical 
models with finite viscosity - The nature of mass evolution is yet 
to be fully resolved. Recent observations of odd— even effects in 
the primary mass and charge yields further complicate the problem- 

Cl. 5] ANGULAR DISTRIBUTIC»C OF THE FISSION FRAGMENTS 

In fission induced by an energetic projectile, the fragments 
exhibit an anisotropy with respect to the incident beam direction. 
Bohr applying the ideas of the collective model to the deforming 
fissioning nucleus, explained the observed fragment angular 
distributions on the basis of his chann&l tfieory .In this theory, 
the angular distribution of the fission fragments are related to 
the quantum levels of the fissioning nucleus atthe saddle point- 
Thus, spectroscopy of highly deformed nuclei comes naturally from 
the study of fission fragment angular distribution measurements. 
Extended to the medium excitation energy region with a statistical 
distribution of levels at the fission saddle point, the measured 
anisotropies yield direct • information on the nuclear moments of 
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inertia and therefore on the shape of the saddle points Underlying 
the Bohr model and its extension to high excitation energies is the 
assumption that the tilting mode attains equilibrium at the fission 
saddle point and does not change during the descent to scission- 
Extensive measurements^^ of fragment anisotropies for a number of 
systems over a wide range of bombarding energies have clearly 
demonstrated the smooth transition from the discrete level regime 
close to the fission barrier to the statistical region, a rapid 
washing out of nuclear shell effects on the level densities with 
excitation energy and shape parameters of fission transition state 
nuclei at high excitation energies consistent with the predictions 
of the liquid drop model- Fragment angular distribution is one of 
the well understood features of the fission process at present and 
it promises to be a senstive probe of angular momentum transfer in 
heavy ion induced reactions- 

A very interesting correlation between the fragment anisotropy 
and the mass has been seen in a few experiments^^- While the 
correlation can be explained in some cases on account of 
multichance fission, the effect appearing to be genuine in the 
others- If so, this demands a reevaluation of our present 
understanding of fragment mass and angular distributions, for, it 
is generally believed that the angular distributions are frozen at 
the saddle point while the mass and charge distributions are 
decided only at the scission point- This may rule out any 
correlation between the two observables- 
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11.61 SECONDARY RADIATIONS FROM FISSION FRAGMENTS 

The configurations of the two fragments touching each other at 
scission determines various distributions in fission- It is 
beleived that both the bulk properties as well as the microscopic 
properties are decided during this stage of motion of the nuclei on 
its way to fission- At the instant of separation, the primary 
fragments acquire considerable energy from the following sources- 
the intrinsic excitation energy of the fragments at scission, the 
deformation energy of the fragments at scission and the coulomb 
excitation of one of the fragments by the other during separation - 
The fragments then deexcite through the emission of y neutrons 
and sometimes a charged particle- By studying the mass, charge and 
kinetic energy of the fragment and the properties of the secondary 
radiation by which fragments deexcite, one can study the nature of 
the collective motion from saddle to scission- Most of the 
excitation energy is released in the form of neutron emission and 
so studies of prompt neutron emission provide information about the 
deformation and excitation energies of the fragments at scission- 
The majority of the neutrons are emitted from fission fragments in 
a very short time <less than 4 x 10 seconds after scission)- The 
number of neutrons as a function of fragment mass for thermal 
neutrons was first observed by Fraser and Milton**- The number of 
neutrons <i>) shows a typical saw tooth structure first increasing 
with mass upto symmetric division, then falling sharply beyond 
symmetry and after that increasing with mass- It has also been 
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observed that with the increase in excitation energy of the 

43 

compounrd nucleus this saw tooth structure gets washed away « This 
behaviour has been interpreted on the basis of shell effects and 
deformation characteristics of fragments- Thus^ neutron 
measurements are expected to yield, apart from the magnitude of the 
excitation energy that is transferred into internal excitation 
energy, also the information on the role played by the shell 
deformation and pairing effects on the deexcitation of the 
fragments- 

Studies of photons were stimulated in the 1950’s by the 
realization that the experimentally observed gamma ray energy from 
fission is greater than the theoietital estimates based on the 
assumption that gamma emission takes place after neutron emission- 
The tota^ energy going towards gamma ray emission is about 8 MeV/ 
fission^^- The average binding energy of the neutrons in these 
neutron rich fission fragments is around 4—5 MeV- Thus, the gamma 
rays taking away 8 MeV is slightly puzzling as one more neutron can 
be emitted- One of the reasons could be that the fission fragments 
carry some angular momentum also and the gamma rays are emitted to 
dissipate this angular momentum- Calculations including angular 
momentum effects predict the value closely but they predict that 
above the neutron binding energy the photon emission approaches 
zero- 

The angular distribution of gamma rays is forward peaked with 
respect to fragment direction- Strutinsky^"* ha» interpreted this in 
terms of aligned scission dynamics resulting in fission fragments. 



23 


even from even— even nucleus in its ground state, possessing 
appreciable angular momentum- The gamma rays from fission 
fragments exhibit an anisotropy of 12 to 15 % and it is dependent 
on the energy of gamma ray- The study of anisotropy in gamma 
emision is also important to investigate if the energy from fission 
channel is being transferred into the giant resonance- The fission 
fragments are known to possess 10— 15fi of angular momentum and 
15—20 MeV excitation energy- At this energy most of the angular 
momentum is taken away by the gamma ray- The number of gamma rays 
that are emitted in fission depends not only on the excitation 
energy of the fragments but also on the nature of the excited 
states of these neutron rich fragments- The multiplicity studies in 
heavy ion reactions are known to provide the spin distribution of 
the emitting nucleus- Although there exists a huge data bank with 
regards to the average number of gamma rays emitted in a variety of 
nuclei, there are only a few investigations on the average width of 
the gamma ray multiplicity distribution- The average width is an 
important parameter in the estimation of the spin distribution of 
the fission fragments- The width in the angular momentum 
distribution and even the relative spin that is imparted to the 
fission fragments arises due to the saddle to scission dynamics- 
Thus, the study of secondary radiations which occur at time scales 
orders of magnitude greater than the saddle to scission dynamics 
throws some light on the nature of dynamics involved. 

A detailed study of the multiplicity distribution is called 
for since the coricepts of single particle degree of freedom, the 
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nuclear- shape degrees of freedom, angular momentum and pairing all 
influence the gamma ray emission process very crucially- All these 
various shell effects are predicted to decrease and then vanish at 
sufficiently high excitation energy- Gamma rays from fission 
fragments with moderate amount of excitation energy and spin 
therefore provide an excellent tool in this endeavour, because a 
rapidly rotating nucleus is quite sensitive to shell effects, 
changing rapidly with angular momentum and deformation. Further, 
all these effects are present in the phenomenon of fission in which 
the density of levels in a highly deformed and often rapidly 
rotating system governs the fate of the nucleus- It is in these 
limits of exciation energy and angular momentum that our 
experimental knowledge is very limited- We discuss more of this in 
chapter 2 wherein we have calculated the spin distribution of the 
fission fragments resulting from scission dynamics- We also present 
the variation of the width and. mean multiplicity as a function of 
charge ratio of the fission fragments- We also present our results 
on the prompt gamma ray anisotropy as a function of gamma ray 
energy. The anisotropy studies are important in the study of the 
competition between the El and the E2 modes of emission- This 
competition has been seen in the deexcitation of the heavy ion 
reaction residues and has prompted us to investigate the nature of 
deexcitation of fission fragments to look for this competit ion- 
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[1.71 CONCLUSI ONS 

The fission phenomenon offer many opportunities for the study 
of a phenomenon occuring in nuclear matter under extreme conditions 
with respect to shape, spin and excitation energy- It is thus clear 
that the fission process probes very nearly all aspects governing 
these lax-ge scale collective motion- the structure of the 
multidimensional deformation potential energy surface, the inertial 
and the dissipative aspects of the collective dynamics, and the 
general nuclear behaviour far off the equilibrium configurations. 
While heavy ion reactions have considerably widened the scope of 
such studies, the simplicity of the fission process has retained 
its elementary character. There still remain several aspects of the 
fission process which continue to defy a full understanding even 
today- Fifty years after its discovery, the fission process 
continues to be unique and challenging- 
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CHAPTER II 

GAMMA EMISSION FROM FISSION FRAGMENTS 

12.11 INTRCMXICTION 

Gamma ray emission by excited nuclei, be it from fission 

fragments or reaction products formed in nuclear collisions or a 

product of radioactive decay, has been a rich source of information 

on nuclear structure and reaction mechanism* In fact, nuclear 

spectroscopy as a well defined branch of nuclear physics grew 

largely as a result of measurements of discrete gamma ray 

transitions close to the ground state* More recently, studies of 

discrete gamma ray transitions in rapidly rotaing nuclei has led to 

1 2 

the discovery of backbending phenomenon ' , in particular, and more 
generally to a detailed information on the structure of nuclei with 
large spin* j 

There is also a continuum component to the gamma ray spectrum 
arisisng from the statistical nature of gamma ray emmision at high 

■ I 

excitation energies of the emitting nuclei- Studies of the j 
continuum component , however, are hampered by the unavoidable 
interference from the compton continuum- Nevertheless, continuum | 
gamma rays have provided valuable information- For example, in j 
fission fragments it is the gamma ray measurements whi ch led to the 
informatioTi on the magnitude of fission fragments spin^'^- It is 
now known that even in the case of spontaneous fission of 
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everf-even nuclei fission fragments carry about 10—15 units of 
angular momentum arising due to scission dynamics- 

In the last few years, there is also an increasing interest in 
the study of continuum gamma rays from nuclei with large spin 
produced in heavy ion reactions- These have now become feasible 
with the development of large volume germanium detectors with 
anti-compton shields which results in appreciable reduction in the 
Compton background gamma rays- The study of the continuum gamma 
rays has opened up new vistas in the study of the mechanism of 
gamma ray emission from hot rapidly rotating nuclei- In light of 
some recent results from the heavy ion laboratories, we have had 
to revise a few of our earlier concepts- It is now conjectured^'^ 
that the El statistical gamma rays and the E2 rotational cascade 
compete with each other, whereas previously it was believed that 
the nucleus reaches the yrast state by the emission of statistical 
gamma rays and then dissipates the angular momentum and rotational 
energy by the emmission of E2 cascade along the yrast line- 

One of the ways to study this possible competition between the 
El and the^ E2 modes of gamma emission is to investigate the nature 
of gamma rays emitted from the fission fragments- This study has a 
distinct advantage as compared to the study of gamma rays from 
reactions either induced by heavy ions or by light ions- This is 
because in case of light ion induced reactions, most of the energy 
brought in goes into heating the system giving it very little 
angular momentum, thereby populating the system at the extreme left 
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Fig. 2.1, A typical yrast plot showing the region populated by 
heavy ion and light ion induced reaction as compared 
to region populated by fission fragments. 
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of the yrast plot- Fig- [2-1] illustrates the region populated by 
the light ion induced reaction on a typical yrast plot- The 
nucleus, in this case, is expected to decay primarily by particle 
emission and once this is inhibited then by the emission of El 
statistical gamma rays- On the other hand, in reactions induced by 
heavy ions, alongwith the excitation energy, a lot of spin is also 
imparted to the nucleus thereby populating it upto the extreme 
right of the yrast plot as shown in the Fig- [2-1]- The nucleus is 
spinning rapidly and a substantial amount of the energy brought in 
is locked up as rotational energy- Though the nucleus is expected 
to decay both by E2 rotational cascades and El statistical gamma 
rays, there is uncertainity associated with the reaction dynamics 
and the nuclear matter is thermodynamically hot- The single 
particle properties and the pairing effect is therefore expected to 
be considerably diluted- 

The fission fragments with spin of 10—15 fi and excitation 
energy of 15—20 MeV are in the middle region of the yrast plot as 
shown in Fig. [2-1]- The competition between the E2 and El gamma 
rays is expected to be enhanced. The role played by close shell 
structure, pairing and deformation in gamma emission is better 
investigated in case of gamma emission by fission fragments- We can 
study nuclei with Z ranging from 35 to 65 in a single experiment 
•Jnder identical conditions- This spans the region of nuclei which 
*ne spherical, deformed and soft- In case of fission fragments, 
another advantage is the gamma spectroscopy of neutron rich nuclei, 
which is beyond the realms of heavy ion accelerators but comes 
naturally in this case- 
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The fission fragments emit gamma rays in the last stage of 
their de— excitation process^- The fission fragments carry a lot of 
excitation energy arising due to scission dynamics and de— excite 
first by neutron emission which takes place within about 10 to 

seconds- Gamma ray emission takes place after that, upto a 
time as late as 10 ^ seconds- There are also delayed gamma rays 
which are mainly due to radioactivity of the neutron rich fission 
products and isometric transitions- These delayed gamma rays 
contain very little information about the fission dynamics and we 
will discuss them no further- 

From the study of angular distributions and the multiplicity 
distributions of prompt gamma rays, it is possible to infer the 
spin of the fragments- The overall picture, now, is one in which 
the fragments in spontaneous and thermal neutron induced fission 
appear to have 10—15 h of angular momentum oriented preferentially 

o ^ 

about 90 to the fission axis - The electrostatic forces betwen the 
fragments at scission would induce such angular momenta in the 
fragments if the fragments are formed in axially asymmetric 
configurations. Although the de— excitation mechanism of the fission 
fragments has been a subject of many investigations in the past, 
information regarding gamma ray emission, such as its dependence on 
the single particle properties, pairing, deformation and spin 
distribution, is still lacking. 

The experimental observables in gamma ray measurements are the 
first moment <N^>, and the higher moments <N^^> of the multiplicity 
distributions- The first moment is the average number of gamma rays 
emitted by the fragments, the second moment is related to the width 
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of the di s tribut i on and the third moment gives the skewness of the 
distribution^^- On each one of these, a differential measurement 
with respect to the mass, charge or kinetic energy of the fragment 
can be made- Alternatively, one can measure these quantities 
averaged over any one of these quantities or over all of them- One 
can also measure the anisotropy in gamma ray emission with respect 
to the fragment direction averaged over the fragment properties or 
differentially as a function of any one of them- 

In this chapter, we present our work which aims to measure the 
average multipicity and the width of the multiplicity distribution 
of gamma rays as a function of the charge ratio of the fission 
fragments- We have also measured the anisotropy in gamma ray 
emission With respect to the fragment direction as a function of 
gamma ray energy- Using the available expressions for the branching 
ratios for neutrons and gamma ray emission from fission fragments, 
we have written a statistical model code based on the Monte Carlo 
technique which traces the history of the de—exci tation of the 
fission fragments which have excitation energy of about 15—20 MeV 
and angular momentum of 10—15 fi» 

This chapter is divided in five sections- In Sec- [2-23$ we 
describe the details for the measurement of the multiplicity and 
the width of the multiplicity distribution of the gamma rays 
emitted from fission fragmentsCin spontaneous fission of ^^^Cf) as 
a function of their charge ratio. In Sec.C2.3], we present the 
anisotropy in prompt gamma rays from fission fragments of ”^Cf 
with respect to the fragment direction of motion as a function of 
gamma ray energy. In Sec- [2-43, we describe the Monte Carlo based 
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statistical model code to simulate the gamma ray and neutron 

emission from the hot rotating fission fragments. Here, we 
investigate the role played by shell structure, excitation energy 
and the angular momentum carried by the fission fragments on the 
mean multiplicity and the width in the gamma ray multiplicity 

distribution in the spontaneous fission of Cf - Lastly, based on 
our theoretical calculations, we have made an attempt to calculate 
the initial spin distribution of the fission fragments- Sec- [2-5] 
contains the conclusions drawn from the above studies- 

[ a* 2 ] MULTI PLI Cl TY MEASUREMENTS 

C2-23<a> Introduction 

There have been a number of measurements in low energy fission 
on the average multiplicity of the gamma rays and the average 

energy taken by them in order to investigate the nature of 

i.1. 12 13 

de— excitation of the fission fragments ' - Pleasanton measured 

the multiplicities of the gamma rays in thermal neutron induced 

fission of U in the time range of 5 ns, 70 ns and 275 ns- His 

value of the multiplicity agree well with the value obtained by 

Verbuisky^*. This implies that most of the gamma rays are emitted 

12 

within 5 ns of the fission process. Mehta etal- measured the 
average multiplicities and the average energy taken by the gamma 
rays in the ternary fission of Cf- The average energy taken by 
the gamma rays is not much different than in the case of binary 
fission- There have been some measurements^^ which are in 
contradiction to the results of Mehta etal -- Increased RT indicates 

■ y ■ 



that the spin of the fragment is more but in ternary fission, due 
to stretched scission configuration, the spin should be less- 
However, Ramamurthy etal-^^'^^ have also found that the FJ^ does not 
show any significant variation in ternary and binary fission- This 
may imply that the average excitation energy and the spin of the 
fragments remain the same irrespective of the emission of light 
charged particles- These results are well documented and the 
average energy taken by the gamma rays is 7-0±0-5 and the average 
number is 9.7±0-4^^. 

It has also been observed that there is a positive correlation 
between the number of neutrons and the average number and the 
energy of the gamma rays emitted by the fission fragments while 
de— exciting. Johansson observed that the gamma ray multiplicity 
as a function of fragment mass shows a prominent saw tooth 
structure which is similar to the structure exhibited by the 
neutrons - The gammma ray multiplicity increases linearly with 
fragment mass till the doubly closed shell at mass 132- It then 
falls sharply and then again increases linearly with the fragment 
mass- In case of neutrons, this saw tooth stucture is destroyed as 
the energy of the fissioning nucleus increases. However, 
experimental information on the variation of multiplicity with the 

excitation energy in the case of gamma rays does not exist- 

20 

Nifenecker obtained a closed form expressions for the correlation 
between the number of neutrons emitted and <i>the average energy 
and (ii)the average number of gamma rays respectively! 

<E Cm,E >> = 0-75 v<m,E > + 2 

K ’ K 

<N^<m,E >> = 1.1 KmyE ) + 1.75 

iC 
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where E is the average energy taken by the gamma rays • 

is the average number of gamma rays emitted* 
u is the average number of neutrons emitted- 
m is the mass of the fission fragments- 

is the kinetic energy of the fission fragments- 
However, all these measurements have been made as a function 

j 

of fragment mass which is obtained by measuring the energy of the 
two correlated fission fragments and then invoking mass and 
momentum conservation principles'^- The mass resolution obtained by 
this technique is of the order of 4 amu- The pairing effect, if 
playing a role in the gamma emission process, is expected to be 
absent- The deformation effects are also expected to be diluted 
becuse of the poor mass resolution- 

There has only been one measurement on the width of the gamma 
ray distribut ion^^'^^^ where the width was measured as a function of 

2!52 

single fragment kinetic energy in the spontaneous fission of Cf- 
The measurements were done both for binary as well as ternary 
fission^^- The results indicate that width in case of ternary 
fission is larger than that in binary fission- This, according to 
the authors, is because the average spin in case of ternary fission 
is small leading to less number of cascade gamma rays- However, in 
measurements as a function of single fragment kinetic energy the 
shell properties are expected to be diluted- The second moment of 
the multiplicity distribution as in heavy ion reactions is expected 
to yield information about the original spin distribution of the 
nucleus • The present experiment was performed with the aim to see 
the effects of pair correlations and deformations on the average 



multiplicity and the width of the multiplicity distribution- The 
charge of the emitting fragment was identified by measuring the 

characteristic K X-ray it emitted^^- By this technique, we could 
get a single charge resolution. 

[2.23<b> Experimental Details 

A californium 252 source with a strength of 5 /j Curie was 

mounted in a vacuum chamber. A solid state surface barrier detector 
was kept at a distance of 3 cms- from the Cf source to detect 
the fission fragments- A high purity germanium <HPGe) detector was 
kept in a direction perpendicular to the direction of fragment 

motion to detect the K X— rays- The distance of the HPGe detector 

was around 4 cms from the source- The Californium source was 
mounted in such a way that it subtended an angle of 45*" to both the 
fission detector and the HPGe detector- Two sodium iodide (Nal> 
scintillation detectors were kept in the plane perpendicular to the 
X— ray direction- One of the sodium iodide detectors was a 2” X 2“ | 

crystal and the other was a large 5" X 5" crystal- The two 

scintillation detectors were kept such that one subtended an angle | 

of 160* to the direction of fragment motion and the other an angle | 

of 200* to the direction of fragment motion- Fig- [2-23 shows a 

schematic arrangement of the experimental set— up- The large sodium | 
iodide detector was suitably shielded with lead bricks to avoid the 
detection of stray neutrons and cosmic ray particles- Before 

starting the experiment, it was ensured that this shielding was 
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For fission 



Fig. 2.2. Geometrical setup for the measurement of the multiplicity 
and width of the multiplicity distribution of prompt gamma 

rays emitted from fission fragments in spontaneous fission 
of 252 q^ . 



effective by recording a K X— ray spectrum in coincidence with 
fission and gamma ray without the Californium source- 

To measure the average number of gamma rays emitted per 
fission, it is sufficient to employ just one gamma ray detector of 
known detection effeciency and of known solid angle factor for a 
fixed number of fission events- To measure the width of the gamma 
ray distribution, it is essential to employ one more gamma ray 
detector- In fact, the number of detectors needed is directly 
related to the moment of distribution one is aiming to measure- 
This is made clear by the set of equations that we present in the 
next paragrap»h- 

The equations given below connect the coincidence rates to the 
multiplicity distribution P 

n 

Ci= ^ Cl - < 1 - Oi <2 .i> 

i = l,2 

c = r Cl - c 1 - O.)" - < 1 - 0 . >"+ < 1 - 0 . -O. >" ]P <2. 2) 

12 1 j 1 j n 

i = l,2 ; j==l,2; i = j 

where O. is the detection effeciency and C. is the coincidence rate 
of the i’th gamma ray detector- includes both the detection 

effeciency as well as the solid angle factor- If we assume that the 
solid angle of gamma ray detectors is small and the probability of 
two gamma rays entering one detector simultaneously is small, then 
for a single fission event the number of gamma rays detected will 
be given by: 

C.=<ri >0. 

1 y ^ 

where <n^> is the average multiplicity- 
The probability that after the first gamma ray detector has 
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detected one of the gamma rays, the seco ^etector will 

detect one of the <ri^“-l> gamma rays is : 

C. . == < 1-1 <n ~1>>0. 0. 

ij T' 2^ 

It is therefore seen that the count rate for single and double 
coincidence are related to <n > and <ri^>« Since Q,s are vne only 

r r ^ 

unknown quantities, to eliminate them it is sufficient to have a 

priori knowledge of n^« (However, since we did not do an absolute 

efficiency calibration of our gamma ray detectors, we have taken 

the average number of gamma rays n^ from the data by Skarsvag - He 

252 

has measured it for spontaneous fission of Cf and according to 
him fT =9-7. In writing the expression for <n > and <n^> we make two 

r r r 

assumptions : 

(i>the two gamma rays are not correlated 

<ii>the detection effeciency of our gamma ray detectors averaged 
over the gamma ray spectrum is nearly constant for single and 
double coincidence measurement s - 

Assumption <i), in light of the fact that only about 5 gamma 
rays in 411 geometry are emitted from one of the fragments, is 
approximately correct- Assumption (ii> is also justified in light 
of the fact that the effeciency curve for the Nal detector for 
gamma ray energy below 5 HeV is almost flat and the number of gamma 
rays above 5 HeV emitted by fission fragments is very smal 1 - 
Figure 2-3 shows the block diagram of the electronic set up 
used for this experiment- Pulses from the HPGe detector i-e- the 
K— Xray pulse was fed to a preamplifier and then to a linear 
amplifier- The pulse from the linear amplifier was fed to the input 
of the linear gate and stretcher- The fission fragment pulse from 




Fig.2-3. Electroriic stit up for the «letermination of the 
multiplicity and width in the mulplicity of gamma rays 
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from fragments in spontaneous fission Cf- 









the solid state surface barrier detector was passed through a 
preamplif ier-amplif ier combination and then to a single channel 
analyzer, which was used to cut off the natural alpha particles- 
The pulses from the two gamma ray detectors were also processed 
similarly- The single channel analyzer in case of gamma ray 
detectors were adjusted to cut off gamma ray energies below 80 keV- 
The pulse from these single channel analyzers were fed to a 
mult i coincidence uni t which had the option of switching off or on 
any one of these pulses- 


C2.23<c) Results 

First the K—Xray spectrum was recorded in a multi channel 
analyzer in coincidence with fission fragments- We will call this 
spectrum the double coincidence spectrum, which is shown in 
Fig- [2-4]- The triple coincidence spectrum was recorded by 
switching on the mult i coincidence unit by the gamma pulse from 
either one of the Nal crystals alongwith the fission fragments- The 
triple coincidence spectrum obtained by the larger sodium iodide 
crystal is shown in Fig-C2-5] Finally, we recorded the K“Xray 
spectrum in coincidence with the pulse from the two gamma ray 
detectors and the fission detector- This we will call the four fold 
coincidence spectrum- The four fold coincidence spectrum is shown 
in Fig. 1 : 2 -5]. 


These spectra were then unfolded to get the intensity due to 

various charges- For this, a fitting program was used which 

required the energy and the relative intensity of the K and X- 

a ft 
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ray lines of the various charges - The program also required the 
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dstector response function* The background was fitted with a second 
order polynomial fit* x minimisation was done to get the 
proper unfolding of the spectra* The area under each of the charges 
gives the contribution of that charge for either the double, triple 
or four fold coincidence spectra* From the measured K X-ray 
intensities due to various charges in the double, triple and. four 
fold coincidence spectra, we have derived the average gamma 
multiplicity <N^> and the width of the gamma multiplicity 
distribution as a function of the charge ratios of the fission 
fragments* The procedure followed is simple but for the sake of 
completeness we will describe it* For a particular charge if N* is 
the average number of fission events that the fission detector has 
detected; N* is the average number of X— rays emitted peT* fission 
and is the average number of gamma rays emitted per fission, 
then the intensity of that charge in the double triple and four 
fold coincidence spectra is given by 
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Fig- C2-7] shows the intensities obtained due to various charge 
ratios <C^, C^- C^) in the double triple and four fold coincidence 
spectra. 

Fig. [2-8] shows the average multiplicity as a function of 
charge ratio of the fission fragments. The multiplicity seems to 
indicate some odd even differences, for example, multiplicity for 
charge 54 is larger than that for 55 and again it is larger for 
charge 56 than for 57. However, since the error bars are large any 
definite statements cannot be made on average multiplicity as a 
function of fragment mass. 

Fig«C2-9] shows the width of the multiplicity distribution as 

I 

a function of the charge ratio of the fission fragments. The width 
clearly shows a dependence on the odd even nature of the fragments. 
The widths corresponding to charge 52 is larger than for charge 53, 
for charge 54 it is again large and then for charge 55 it is small. 
As we pass over to the region of deformed nuclei <Z= 57— 59>, the 

odd even character seems to vanish. Slightly lower values of the 
width observed and the vanishing of the odd even structures can be 
attributed to the fact that the nucleus in this region is a rapidly 
spinnig object. We can deduce from our experimental observations 
that the magnitude of the rotation is such as to destroy the 
pairing correlat ions . Secondly , a substantial amount of energy is 
now stored as the rotational energy of the fragment- The fragment 
dissipates this rotational energy and angular momentum by the 
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Bar chart showing the intensity of 
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Fig. 2.8 Variation of multiplicity of the prompt gamma 
rays as a function of charge ratio. 
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Fig. 2.9. Variation of width of the multiplicity distribution 
as a function of charge ratio. CEMTRsl j 
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emissiori of E2 rotational cascades- These are transitions from 
fixed energy states with very - little statistical variations 
associated with them- The contribution to the width comes from the 
'statistical nature of emission of neutrons and the dipole 

transitions- The sp«herical nuclei are expected to de—excite by the 
emission of neutrons and then by statistical gamma rays- Being 

rotationally degenerate systems, the emission of E2 rotational 
cascades is inhibited- Thus the width of multiplicity distribution 
for spherical nuclei is expected to be more than that for a 

deformed nuclei- This is supported by our experiment also- 

If we assume that the nuclear temperature at the scission 
point is 0-5 to 1 Mev and the moment of inertia of the fission 
fragment is between 100-200 MeV- Then the value of the experimental 

.’N 

width obtained is in general agreement with the value of the width 
obtained by Moretto^"* assuming that the energy is partitioned 
equally amongst the various angular momentum bearing modes- The 
angular momentum bearing modes are the <i> bending mode <ii> 
tilting mode <iii> twisting mode and (iv> wriggling mode- The 



of this chapter after we have presented our results for the 


statist! cal model calculations- 



54 


£2.31 ANISOTROPY STUDIES 

J 

C2-33<a> Introduction 

I The angular distribution of the y rays emitted from the 

! I 

fission fragments with respect to the fragment direction can give 
information about quantities such as the multipolarity of radiation 
and the alignment of the angular momentum of the fission fragments. 

An angular correlation between the fission fragments and the prompt 
y rays has been observed, with the emission in the direction of the 
fragments more probable by about 15X than emission perpendicular to 
the direction of the f ragments^'^. There are two reasons for the 

anisotropy in the fission gamma ray angular distribution* One is 

\ 

\ 

due to the Doppler effect i.e. an isotropic distribution of f rays 
in the frame of reference of the fragments will give rise to a 
distribution in the laboratory system due to Doppler effect given 
bys 

W<©> = W Cl+<2+r>v/c cos©] 

where v is the velocity of the fragment, c is the velocity of 
light, © is the angle with respect to fragment direction, and r a 
small correction factor accounting for the change in energy of the 
y rays. The Doppler anisotropy disappears if the experiment does 
not distinguish between the fragments or if the fragments are 
stopped before they emit y rays- 

The other cause of the anisotropy is of our interest. It is a 

* 

consequence of a preferential orientation of the fragment spin with 
their direction of flight- The nature of anisotropy is a pointer to 
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the fact that y emission is predominantly of quadrupole nature- Any 
angular momentum present in the original fissioning system appears 
later in the final fragments <but with the axis of rotation tending 
to lie in the direction of fragment motion) and as relative angular 
momentum between the fragments- The initial angular momentum, for 
low values of the projection of angular momentum on symmetry axis 
<which is expected for spontaneous and thermal neutron induced 
fission) is expected to appear as orbital angular momentum between 
the fragments. It has also been found that the magnitude of the 

i 

anisotropy is irisenisi tive to the initial compound nucleus spin* 
This along with the non appearance of the initial compound nucleus 
spin in the orbital angular momentum implies that the origin of the 
anisotropy with respect to fragment direction lies in the descent 
of the nucleus from the saddle to scission, a fact that can be 
utilized to study the dynamics of the fission process- 

25 

According to Strutinsky , axially asymmetric scission, gives 
rise to a transverse component of coulomb repulsion between the 
nascent fragments, which will cause them to rotate after fission in 
opposite directions- He has also derived a simple expression for 
the anisotropy of the radiation- Assuming that the de— excitation 
proceeds statistically, the angular distribution is given by 

W<e> « 1+K( j/ IT> sin e 

where K = + 1/8 for dipole radiation 

= — 3/8 for quadrupole radiation 

Here,j is the angular momentum of the fission fragments 
I is the moment of inertia of the fission fragments 



56 


T is the nuclear temperature* 

There is always an admixture of dipole and quadrupole 
radiations* The assumption of pure quadrupole radiation leads to an 
under estimation of spin* The assumption that the relative 
probabilities of decay to states of higher or lower spins than the 
emitting' state are dependent on a statistical level density 
weighting factor is also not quite correct for the last steps of 
the cascade where there are no final states of higher spin 
available* 

Wilhelmy^^, by measuring the angular distributions of several 
2' ' > Q transitions in the ground state band of even fission 
isotopes, had established that there is a preferential alignment of 
fragment spin with respect to their direction of flight- They found 
a preferential emission of E2 radiations along the direction of the 
fragment motion with anisotropies ranging between 8*3X to 33-4.X* 
Hoffman*^ had done an early analysis of gross angular distribution 
of the fission y rays* Wilhelmy’s result agrees with that of 
Hoffman that the initial spins of the fragment are preferentially 
aligned perpendicular to the fragment’s path* 

Ivanov has studied the fission y rays as a function of 
fragment kinetic energies, mass ratios and also as a function of y 

23S» 

ray energy inc slow neutron induced fission of U- They found 
that the y anisotropy as a function of total kinetic energy shows a 
definite increase for all mass ratios- But y anisotropy averaged 
over total kinetic energy is rather insensitive to mass ratios* 
Ambruster has studied anisotropy of y rays, emitted between 10 
and 100 ps after the fission process, as a function of fragment 
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mass* Their results are in conf irmation with that of Ivanov- They 
find that the anisotropy is approximately 20 % for the heavy 
fragment mass whereas it is only 14% for the light fragment mass- 

Skarsvag’s work is one of the few available in literature 
where differential measurement of anisotropy has been done as a 
function of y ray energy. They find that total y rays emitted 
within 12 ns after fission with energies greater than 0.114 MeV is 
9-7 + 0-3 per fission and the total y ray energy released is 

<7-0j;p.3> MeV per fission. They have measured y anisotropy as a 

function of y energy from 0.114 MeV to 3.5 MeV. They find that the 
dipole and the quadrupole components are about equally strong at 
high y ray energies, the dipole moment predominates at low energies 
and the quadrupole component at intermediate energies. 

i 

In the present work, we have measured y anisotropy as a 

function of y ray energy from 50 KeV to 5 MeV. 

C2.3]<b) Experimental Details 

252 

A spontaneously fissile source Cf of "Sfj Curie strength was 
kept in a vacuum chamber. Two solid-state surface barrier detectors 
were mounted in the chamber such that both were at 90* with respect 
to each other and each of them saw the source symmetrical and at 

nearly equal solid angles. A NaI<Tl> scintillation detector was 
used to detect the y rays and was kept at a distance of 60 cms from 
the source in line with one of the fission detectors, henceforth 

called as the F detector, and at an angle of 90* to the other 

o 

fission detector, henceforth called the F detector. Thus, the 

PO 

coincidence between the F^ detector and the Nal detector signalled 
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Scintillation datactors 



Fig. 2.10. Geometry of the experimental arrangement for the study of 
prompt gamma ray anisotropy as a function of gamma ray 
, energy in spontaneous fission , of • 
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y ray emitted in the direction of the fission fragment, while the 

coincidence between the F detector and the Nal detector signalled 

oo 

a y ray emitted perpendicular to the direction of fragment motion- 

The fission detectors were ORTEC F series heavy ions 

detect ors« Each detector had an area of 300 with a depletion 

»l u 

depth of 60 p- The sodium iodide detector was of 5 x 5 size with 
a resolution of 10% at 600 keV* Since the neutrons emitted from the 
fragments are strongly correlated to the direction of the fragment 
motion and thereby show a large anisotropy, it was necessary to 
eliminate them by the time of flight technique- Since we had a very 
large Nal crystal, it had to be properly sheilded from stray 
neutrons and cosmic rays- The shielding was tested with a run 
without the source and the background contribution was found to be 
negligible- 

fast coincidence was taken between either of the fission 
pulses and the y ray pulse- To do this, the timing pulse from the 
preamplifiers of the fission detectors were passed through a fast 
timing filter amplifier, the output of which was fed to a constant 
fraction di s criminator which was operated in the external delay 
mode- Walk adjustments were made to optimize the timing 
resolution. The pulse from the dynode of the photomultiplier 
coupled to Nal detector was processed in a similar fashion- The 
outputs of the two constant fraction discriminators being fed by 
the F and F pulses were ORED and were used as a stop pulse for 
the Time to Amplitude converter, the start pulse being given by 

I 

the y ray pulse- • 

The energy pulse from the two fission detectors was passed 
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through a research amplifier which could give a — 2\/ maximuii) signal 
(a requirement of the ADCs). The energy pulse from the y ray 
detector was passed through a Spectroscopy Amplifier. Suitable 
delays were introduced to ensure coincidence between the analog 
pulses and the gate. The electronic block diagram used for the 
experiment is shown in Fig. [2.11]- The fission counts in the two 
fission detectors were monitored after every 12 hours. This was 
done to normalize the data such that both the detectors detected 
equal number of fission events. 

A delay of 50 ns was introduced in the timing channel of F 
detector to seperate the y pulses due to F and F - This resulted 
in a clean separation between the true peaks due to gamma emission 
along the direction of fragment motion and perpendicular to the 
direction of fragment motion. A timing resolution of 2-5 ns was 
obtained. The TAC peaks obtained are shown in Fig- [2. 12]- The 
figure shows the TAC peak due to the F and the F Channels. Also 
shown in the figure are the neutron humps in the direction of 
fragment motion and perpendicular to the direction of fragment 
motion. 

The digital pulse from TAC which was +5 W was fed to a gate 
and delay generator to get a —2 \/, a requirement of the CAMAC A to 
D converters being used- This —2 V logic pulse formed the master 
pulse and gated the camac ADC’s Four parameters, namely the two 
analog pulse of the fission detectors, the cathode pulse from 
photomul t ipl ier carrying information about y ray energy and analog 
pulse from the TAC. These were recorded in the list mode on a 
magnetic tape on a camac based data acquisition systeiii. The 
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1 Electroriic set-up for the determinat ion of anisotropy of 
prompt gamma ray w-r.t fragment motion as a function of 
gamma ray energy- 
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external delay in the F timing channel resulted in clear 

fi >0 

separation of the two true TAG c*eaks. The y spectrum due to F and 
due to were separated during the off line analysis. Above 180 

keV y energy cut off, a clear separation between the true y peak 
and the neutron hump was achieved- 

[2.3] <c> Results 

From the TAG peaks due to gamma rays in the direction of 
fragment motion and in the direction perpendicular to the fragment 
motion with suitable normalization (as defined in the earlier 
section), we could obtain the anisotropy in the gamma emission as a 
function of gamma energy- The nature of our curve is very similar 
to that obtained by Skarsvag - Fig- [2-13] shows the anisotropy as 
a function of the gamma ray energy. We find that in the energy 
region between 200 keV and 700 keV, anisotropy shows a strong 
correlation with the fragment 'direction being a maximum of 
1 •2610-01 at 520 keV- This is a typical yrast bump seen in the 
gamma rays from heavy ion reaction residues too- The value obtained 
by Skarsvag in the energy interval 0-491 to 0-566 is 1-271 0.04. 

The similarity of the curve obtained by us and that by Skarsvag 
extends further also- We have made measurements beyond that done by 
Skarsvag and .we find that the nature of the curve does not change - 
The anisotropy always remains more than one which implies that 
there are more number of E2 transitions than El- This is true at 
gamma ray energies even beyond 2-5 MeV- 

For energies greater than 0-7 MeV, one does not expect any 
rotational E2 transitions- We expect the nucleus to decay by 



Fig. 2.13. Anisotropy in / emission with respect to fragment direction of motion as 
function of /-ray energy. 




neutr^ons or* El trans i t i ons « Since the fission fragments carry a 
moderate amount of excitation energy and spin it is difficult to 
conjecture the emission of E2 vibrational t ran s i t i ons - However , if 
the transitions were predominantly dipole in nature, then by 

2!5 

Strut insky ’ s pres cr ipt i on , anisotropy should have been less 
than one- One conjecture is that these transitions contributing 
towards increase of anisotropy are either E2 vibrational lines in 
the even-even fragments or they are E2 statistical gamma rays- The 
anisotropy shows a similar behaviour in case of heavy ion 

30 3 ± 

data * - The nature of these transitions, there also is not 
known- The angular distributions studies alone cannot distinguish 
the various possibilities and some more experimentation is needed 
to establish the multipolarity of these transi tions - 


[2.41 STATISTICAL MODEL CALCULATIONS 

The nucleus, given a small excitation energy, being a complex 
many body system can exist in many different states- This is 
utilized by the statistical models of nuclear structure and nuclear 
reactions- The fact that the density of quantum mechanical states 
increases rapidly with the excitation energy has not only made the 
statistical models appropriate but also essential for the 
comprehension and prediction of many nuclear phenomena- 

One can use statistical model for calculating the average 
properties of a nucleus and also for the decay of fully 
equilibrated system such as the compoud nucleus- The essential 
assumption in all statistical models is that all possibilities for 



the decay are intrinsically equally likely- Physical process is 
governed by factors such as the density of final states and the 
barrier penetration factors- The probability for a particular decay 
to occur is thus inversely proportional to the total number of 
possilble decays- 

In heavy ion induced reactions, because of high excitation 
energy and spin, the reaction prod»jcts have many decay 
possibilities- The decay possibilities may to some extent depend on 

33 

the entrance channel parameters which are also many - This has 
resulted in a number of computer codes to simulate the statistical 
decay of these hot rapidly spinning reaction products depending on 
the type of information the user wants to extract- These model can 
broadly be classified in two categoriess 

(i) the single step processes <ii>the multistep processes 
In the multistep processes, gamma emission, • particle 
evaporation and fission are included explicitly- The popular ones 

3-4 3!5 3<5 37 38 3S> 

are i4ltce ' » Cascade ' Julian , Roulotte etc- Cfor more 

details see Ref -[33] )- 

All these models have been developed for heavy ion reactions 
where a nucleus is given a lot of excitation energy and angular 
mofnentum- In case of low energy fission, the fission fragments are 
populated with a moderate amount of excitation energy and spin- 
Thus, the decay of these fission fragments is slightly different 
from a nucleus prodouced in heavy ion reactions- To analyze our 
multiplicity data, we have therefore written a statistical model 
code to trace the history of decay of the fission fragments having 
15—26 MeV of excitatioTi energy and around 10-15 units of angular 



momentum- 


Our primary aim in developing the code was to estimate the 
spin distribution of the fission fragments- The hot rotating 
fission fragments first cool by the emission of neutrons and then 
they further cool by the emission of the El statistical gamma rays 
till they reach the yrast lj.ne- After reaching the yrast state the 
fission fragments dissipate angular momentum by the emission of E2 
rotational cascades - 

The fission fragments are expected to de—excite by emitting 
(i> neutrons 

(ii) El statistical gamma rays 
<iii>E2 statistical gamma rays 
<iv) E2 rotational cascades 

Uhile the neutrons and the El and E2 statistical gamma rays 
cool the nucleus, most of the angular momentum is dissipated by the 
E2 collective rotational lines- El statistical gamma rays can also 
take away one unit of angular momentum, but the most probable 
transitions are the ones which leave the nucleus in the same 
angular moihentum state as the initial one- While E2 statistical 
gamma rays can take away two units of angular mometum, their 
emission probability is much smaller. Therefore, most of the 
angular momentum is dissipated by the E2 rotational transitions. 

Monte Carlo procedure, is used to calculate the number and 
energy of the gamma rays and the neutrons emitted by a fragment 
given some excitation energy and spin- Following steps illustrate 
the method briefly! % 

Step 1 ! It calculates the partial decay width for the possible 



decays using the following expr-essions s 

E-B 

r = <D/2n)<:2mR^g)X " £p<E-B -£:>d£ 

r» O n 

M I ^ - 

where D is the level spacing. 

m is the mass of a nucleon- 
R is the radius of the fragment- 
g is the intrinsic neutron degeneracy- 
is the binding energy of the neutron- 
\ is the multipolarity of the radiation- 

is a constant, dependant on the multipolarity of the 
radiation. 

p is the level density parameter; a function of the spin 
and the excitation energy of the fragment- 
The level density parameter used was that of Ignatuk and 
Smirenkin**^ and is given bys 

p(E ) = a<l-/<x)A /E > 

X go X 

where /<k> = [1 — e ^^x] 

a is the asymptotic value of the level density paramater 
and is given by s . 
a = ACoi + 

ot = 0-0730 f /9 = 0-115 

A is the ground state shell correction 
st» 

r = 0-054MeV“‘ 

The spin dependence In the level density parameter was introduced 
by reducing the excitation energy by the amount of energy stored 



as rotational energy. 

Step 2" The partial decay widths for various particles are 
calculated which determine the weight factor for the various types 
of decays. 

S'fcgp 3: Once the type of decay is decided, it determines the energy 
of the transition and reduces the excitation energy of the fragment 
by the amount of energy of the emitted transition- In case of 
neutron emission, the excitation energy is further reduced by the 
bindig energy of that neutron- 

Step 4: The procedure is repeated with reduced excitation energy 
(and mass if a neutron is emitted). 

The iterations are continued till the difference in the 
excitation energy of the fragment and the rotational energy is less 
than 0-5 MeV- The iterations are stopped at this difference because 
below this energy the decay may not be statistical due to scarcity 
of energy levels- In fact, placing a cutoff less than this distorts 
the gamma energy spectrum at low energy. The probability for low 
energy transition becomes very high, giving a spurious peak in the 
gamma ray energy spectrum. 

For calculating the average multiplicity and width in the 
multiplicity distribution, the history of a particular fission 
fragment, having a fixed excitation and spin, was traced a hundred 
times. The same was done for the complimentary fragment having the 
same value of spin and excitation energy. The calculation for the 
multiplicity and width of the multiplicity as a function of charge 
ratios was done by simply adding the average due to each of 
them, while the widths were added in quadrature. 




Charge 

Fig. 2.14. Variation of multiplicity with the charge of the fission 
fragments for two values of Ags (ground state shell 
correction in level density parameter). 



The probability for the emission of E2 statistical gamina rays 
i^as found to be rather small- In a typical case of thousand 
iterations I the frequency of emission of an E2 statistical gamma 
ray was barely 3 or 4- Hence, we did not find it necessary to 
calculate the probability of emission of E2 statistical gamma rays 
and they were removed from the calculations. 

To investigate the role played by the, shell structure in gamma 
emission process, we performed the calculations for the variation 
of <N^> as a function of charge ratio for two values of the shell 
correction parameter A <A=0 and A=4) . Fig.C2.12] shows the average 
multiplicity as a function of charge ratio for the two values of A. 
It shows that the shell correction parameter does not alter the 
average number of gamma rays appreciably. 

The calculations for the average multiplicity and width were 
done .for a wide range of excitation energies and spins. Table- £2. 11 
shows the average multiplicity for excitation energy equal to 20, 
15 and 10 MeV. The calculations were repeated for three different 
values of spin i-e. 15fi, lOfi and 5fi. Table[2.11 shows the average 
multiplicity for the above mentioned values of excitation energy 
and spin for fission fragment pairs having a charge of 56 and 42. 
We find that as the spin of the nucleus is reduced, the 
multiplicity increases. This is expected as the excitation energy 
locked up in rotation is now available for the emission of one more 
gamma ray. The calculation estimates a value of average 
multiplicity which is in close agreement with our experimental 


result 



Table C2- 13 Showing the variation of the multiplicity with the 
excitation energy and spin for the fragment pair 
having charge,«Z<H)=56 & Z<L>=42 


EXCITATION 

ENERGY 

SPIN 

A=4 

AVERAGE 

^MULTIPLICITY 

A=o 

AVERAGE 

MULTIPLICITY 


15 

9. 

9 . P£> 

20 

lO 

4. 44 

5. OO 

5 

3. 58 

3.75 


i 5 

8. 22 

9 . 84 

15 

lO 

a. 48 

8 . PP 


5 

2. 40 

2 . PO 


15 

2. 74 

2.75 

10 

lO 

2, 75 

4. 9 8 


5 

4. 57 

5. 08 


Ue also calculated the width of the fflultiplicity distribution 
as a function of excitation energy for different values of spin as 
for the calculation of average multiplicity* TableC2>23 shows the 
average width as a function of excitation energy again for three 
values of excitation energy i>e-20, 15 and 10 MeV and three values 
of spins i>e> 15fi, lOfi and 5fi> 


TableC2«23 Showing the variation of width of the multiplicity 
distribution with excitation energy and spin for 
pairs having charge Z<H>=56 4 Z<L>=42 


EXCITATION 

ENERGY 

SPIN 

A::^4 

AVERAGE 

WIDTH 

A=o 

AVERAGE 

WIDTH 


15 

1. 19 

1.18 

20 

lO 

1.18 

1. 34 

5 

1.71 

1. 80 


15 

1 . 20 

1.28 

15 

lO 

2 . 40 

2.5 5 


9 

1 . 80 

1 . pa 


15 

0.85 

0 . 8 2 

iO 

lO 

O . P4 

O . P 5 

- 

5 

O. P4 

1.0 2 













The point worth attention is that the average width is very 
much underestimated by the theory as compared to the experiment- 
The difference is not wiped out even for a wide range of excitation 
energy and spin- 

% 

Table C2-33 gives the calculated values of the <N > and < 0 * > 

Y n 

for the various charges- The table reported here is for a spin of 
each fragment being 10 and having an excitation energy of 15 MeV- 


TableC2-3] Showing the multiplicity and the width of the gamma 
ray d istribut ion- 

A = 4 


GHAHOC 

AVCRAOE 

MULTIPLICITY 

WIDTH OF THE 

MULTIPLICITY 


2. 32 

1 . 75 

sa 

2. IP 

1.7 5 

53 

2.37 

1. 55 

S4 ‘ 

2. 37 

1. 54 

55 

2. 75 

1 . 05 

5<S 

2. 57 

1.14 

57 

2. 5i 

1. 01 

53 ' 

2. 70 

1. 18 

5P 

2. 04 

1. 54 

50 

i . pa 

1. 05 


A * 0 


CHAIlOE 

AVERAGE 

MULTIPLICITY 

WIDTH OF THE 

MULTIPLICITY 

51 

2. 54 

1 . pa 

52 

2. 25 

1. 88 

58 

2. 34 

1. 8P 

54 

2.7 3 

1. PI 

55 

3 . 04 

1.22 

55 

2. ap 

1. OP 

57 

3 . 00 

1.08 

58 

3. 13 

1. 25 

5P 

2. 30 

1. 58 

50 

2. 72 

1.0 8 


The salient features of the results of the statistical model 


calculations can be summarised as • 

<i> The probability for the emission of E2 statistical gamma rays 



was found to be rather small. 

Whereas the calculations do reproduce the average 
multiplicity, the width is very much underestimated by the 
theory. 

I 

{iii>We find that the average multiplicity is not a very crucial 
I function of the shell correction parameter A in the level 
density. 

<iv) We find that the average width does not change significantly 
for a wide range of excitation energies and spins. 

[2.51 DISCUSSION AND CONCLUSIONS 

In the light of our calculations showing that the probability 
for emission of E2 statistical gamma rays with tha available 
expression is very small, we are unable to explain the anisotropy 
exhibited by the gamma rays at higher gamma ray energies to be more 
then one- At high gamma ray energies', our experimental results on 
anisotropy indicate that there are more number of E2 transitions 
than El transitions. The possibility of emission of E2 rotational 
cascades of energies greater than 2-5 MeV seems very remote- The 
enhancement of the anisotropy could come from the emission of E2 
gamma transitions, vibrational in nature- Before any conclusive 
statements can be made, one will have to study the multipolarity of 
the these radiations of higher energy. If the nature of these 
radiations is E2 vibrational in nature then one would like to 
investigate as to why theoretically their emission is hindered. 

We have also found that the width in the gamma ray 



fliul i: ip 1 i c i 'fcy distribution is very much underestim-ated by the 
theoretical calculations. The width in the gamma ray multiplicity 
distribution is primarily because of two reasons* 

CD Statistical nature of emission^ of neutrons and gamma rays of 
El nature- Cii) Fission fragments originally having some spin 

I 

distribution. 

If the fission fragments originally have some spin 
distribution, then, depending on the initial spin of the fragments, 
the number of E2 rotational gamma rays emitted to dissipate angular 
momentum will be different. This will then again contribute to the 
width of the multiplicity distribution. 

The width in the gamma ray multiplicity distribution due to 
the initial spin distribution can calculated if we subtract in 
quadrature the contribution to the width due to statistical nature 
of neutrons and El gamma rays from the experimentally measured 
widths. Since we have calculated the width associated with the 
statistical nature of emission of El transitions and neutrons we 
present in the Table(2.4) the width due to the initial spin 
distribution of the fission fragments. _ 



Table(2-4-J The width due to the initial spin distribution of the 
fission fragments- 


GHAROE 

VZDTH IN THE MULT I PL. IC I TV 

DUE TO SPIN DISTRIBUTION 

52 

4. D ±0.7 

53 

3.8 ±0.5 

54 

5.0 ± 0.5 

55 

3.8 ± 0.3 

55 

5.7, ± 0.5 

57 

3.5 ± 0.3 

58 

2.8 ± 0.5 

55 

2.8 ± 0.7 

50 

5. a ±1,0 


The calculated width, as we see from the TableC2-4], is also 
very much different for the spherical fragments from the measured 
experimentally measured widths- The origin of this difference is 
not understood- In our calculations we have assumed that even a 
spherical nucleus dissipates angular momentum by the emission of | 
rotational cascades after reaching the yrast state- There are two 
points we would like to make regarding this assumption* 

<i> The Yrast state is not well defined for a spherical system- 
So the assumption that it dacays along the yrast line is not I 
rigorous- 

<ii) These are rotationally degenrate systems and hence decay by 
E2 rotational cascades should be inhibited- 

However, in absence of any viable theory this was the only 
alternative available at our disposal- In fact, one uses this kind 
of assumption in the statistical de-excitation of the reaction 
products of the heavy ion reactions even if they are spherical- One 
possibility is that they dissipate angular momentum by the emission 
of E2 gamma rays vibrational in nature- This conjecture would 




support our anisotropy results also- 

Our analysis has revealed that the theory in its present form 
is unable to answer the following questions! 

(i) How do spherical and. defopmed nuclei differ in their 
de— excitation mechanism- In particular, by what mechanism 
does a spherical nucleus dissipate angular momentum- 
<ii> How is the angular momentum partitioned between the aligned 
single particles degrees of freedom and collective rotations- 
A theory taking the above mentioned facts into account can 
only explain the mechanism of de-excitation of a hot spinning 
nucleus- 
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CHAPTER III 


THE VARIATION OF THE TOTAL KINETIC ENERGY WITH EXCITATION 
ENERGY IN THE REACTION 

[3.11 INTRODUCTION 

Fission of the nucleus results from a "run axixxy" large 
amplitude collective motion. Even in low energy fission, though 
large energy release is involved, its observables such as fragment 
mass yields and their energies are known to be very strongly 
influenced by shell and pairing effects*. The mass distributions are 
known to be asymmetric in the case of actinide nuclei. 

It is known for a long time that the fragment shell structure 
at scission plays a dominant role in deciding the nature of mass | 
asymmetry. However, attempts to quantitatively calculate the mass I 
and energy distribution on the basis of scission configurations ; 
have not been successful. This is ascribed to the difficulties in | 
the description of the dynamics of the fission process between the 
saddle and the scission point- 

The mass distributions are one of the most important 

charateristi qs of the fission process- Although there is a large 

2 

amount of experimental data on the mass distribution and their 
dependence on excitation energy of the fissioning system for a 
variety of nuclei, a suitable theory to explain the observed! 
characteristics is yet to come. It is known that the actinide | 



nuclei fission asymmetrically leading to two humps in the mass 
distribution- It is also known that the peak to valley ratio in the 
mass distribution is a decreasing function of the excitation energy 
of the fissioning system and at» sufficiently high excitation 
jinergy, the dominant mode of fission becomes symmetric- This can be 
qualitatively explained as due to vanishing of the shell effects- 
The initial excitation energy of the fissioning system, therefore 
plays a crucial role in determining the mass distributions- The 
determination of the mass distribution of the fission fragments as 
a function of excitation energy is important on two counts* 

<i> It tells us as to what role the shell effects play in the 
determination of the mass distribution. 

Cii> While calculating the total kinetic energy as a function of 
excitation energy, it helps us to eliminate the effect of 
excitation energy on the mass distribution by studying the 
total kinetic energy variation as a function of excitation 
energy for mass separated fragments- 

During the last few years, the studies of the dynamical 
aspects of the fission process have rejuvenated, so to say, 
receiving tremendous stimulus from the data from heavy ion reaction 
studies- The discovery of the other modes of large scale collective 
motion, their relation to fission and a unified description is one 
of the arduous tasks before the nuclear physicists today- 

The dynamical aspects play an important role in the last phase 
of the fission process during which the system crosses the fission 
barrier and descends from the saddle point to the scission point- 



present in literature to explain the 


The theoretical models^, 
descent of the nucleus from the saddle to scission, have neither 
been explored in detail nor are sufficiently free from curve 
fitting to be generally accepted^ The models can be broadly 
classified in two groups, depending on whether the motion from the 
paddle to scission is adiabatic with respect to particle degrees of 
freedom or not. 

If the collective motion towards scission is sufficiently 
slow, so that single particle degrees of freedom can easily 
readjust to each new deformation as the distortion proceeds, then 
viscosity plays a significant role in the descent of the nucleus 
from the saddle to scission. In this case, the statistical 
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approximation ' may be valid. Here, it is assumed that the nucleus 
is in statistical equilibrium at all points in the deformation 
space* In that case decrease in potential energy from the saddle to 
scission appears, primarily, as the deformation energy of the 
fragments. On the other hand, if the motion from the saddle to 
scission is fast, then there is no transfer of energy . from 
collective motion to single quasi ' particle excitation- In this 
case, the adiabatic models**"® are valid and viscosity is assumed to 
play an insignificant role in the descent* The difference in 
potential energy from the saddle to scission appears as the 
relative kinetic energy of the system- 

J 

Particle excitations, when dealt with macroscopi cal ly , appear 
in the form of nuclear friction, which is a natural extention of 
the concept of collective coordinates, inertia and potential along 
the lines of classical friction. Hill and Wheeler**, looking for a 



possible explanation of mass asymmetry, in 1953 , suggested a 

microscopic derivation of viscous forces acting during the descent 

from the saddle to scission* In their theory, level slippages 

caused by some residual interaction^ give rise to time dependent 

single particle occupation numbers which, in turn, prevent the 
1 . . . ■ . 

nucleus from fissioning in its lowest state- Their conjecture that 
viscosity may entail asymmetry in the actinide region in low energy 
fission was, however, belied. 

Wilet^*^ gave a full treatment of the interaction betwen 
collective and intrinsic degrees of freedom by time dependent 
Schrtidinger equation- Unfortunately, at that time his idea did not 
receive full attention and it was only by the impact of data from 
heavy ion experiments which brought the idea of deep inelastic 
scattering, fast fission etc- to the forefront that the idea of 
nuclear friction or viscosity received a fresh breath of air- 

Deep inelastic collisions**'** exhibit features very similar 
to those displayed by fission fragments. Here, when a heavy nucleus 
with energy above coulomb barrier and having a large degree of mass 
asymmetry in the target projectile combination, collides with an 
impact parameter b > R, (where R = *^i'*’*^2* *^2 radii 

of the target and the projectile nuclei) with another heavy 
nucleus, it gets scattered inelatically- This scattereing is 
unlike the usual inelastic scattering known in light ion induced 
reaction. ‘ In the exit channel, there is a lot of dissipation of the 
centre of mass energy. The final projectile energy is continuously 
distributed between the bombarding energy of the projectile and the 
fully relaxed value, very nearly equal to the coulomb energy 





between the target and the projectile- The masses of the two 
reaction products are close to the masses of the target and the 
projectile- A naive picture is that, as the nuclear densities of 
jthe two nuclei overlap, damping takes place dissipating the entire 
kinetic energy above the coulomb barrier energy- The two nuclei 
then spin together, when the exchange of a few nucleons becomes 
possible, and after about half a revolution, they separate under 
the coulomb repulsion with very little relaxation in the mass 
asymmetry degree of freedom- The fully relaxed total kinetic energy 
of th& targ&t and th^ proj«c til» lik.& fragments is not given by the 
two spherical nuclei in contact but by those of deformed spheroids 
as in the case of fission fragments- The average time scale for the 
process is around 10 seconds- The energy loss is strongly 
correlated with the scattering angle*®’*'*. There were some early 
attempts to explain the viscous effect*®. Wi lezynski*'* suggested an 
interpretation of, then newly discovered, deep inelastic collisions 
in terms of negative angle' scattering and frictional forces- The 
deep inelastic collisions is a subject in itself and we would not 
refer to it anymore- Our objective was only to show that there are 
other examples of large scale nuclear motion which are viscous and 
for a more detailed study, the reader is referred to Ref- [163 and 
references therein. 

In the case of fission fragments, viscous forces are expected 
to affect the fragment kinetic energies- The fission fragment 
kinetic energy has been measured' from compound nucleus ****Rh 
XNamboodiri*** > upto ®^®110 [Sikkeland*^ , Borderie*® >- If plotted 
against Z*/A*^®, they follow systematically a straight line- This 



behaviour is expected for the lighter compound nucleus upto 135 

2 1./3 

or Z /A 2 s£ 650 because the saddle points of nuclei with fissility 
smaller than approximately x ^ 2^/50A = 0-5 almost coincide with 
their scission points and hence,, the total kinetic energy of 
symmetric fission fragments is well approximated by the mutual 
coulomb interaction energy of two spheres in contact, however, the 
trend seems to be valid even beyond. Fig- [3-1], taken from 
Ref. [19], shows the experimental points of total kinetic energy as 

J 

2 ±/3 

a function of Z /A <a property of the fissioning system)- The 
figure also shows the results of theoretical calculations ' • The 

calculations have been made both by assuming that the system is 
viscous and also without putting any viscosity in the calculations- 
The nature of viscosity assumed for the calculations was of I 
ordinar-y hydrodynaniical Atlnd- While looking at these calculations 
one might say that the nuclear matter is only slightly viscous, the i 
heavy ion data is a pointer to the fact that the viscosity is ; 
large- Thus we can conclude say that the analysis is not unique- In 
fact, Prakash et-al^- have shown f;hat as the fragments separate at 
scission there is considerable dissipation of relative kinetic 
energy in the initial stages due to the proximity friction and a | 
free separation in the coulomb field starts only after the surface I 
separation distance exceeds about 2 fra- This gives rise to a > 
kinetic energy value smaller than the initial energy at scission- 
With the inclusion of this effect they could explain Viola 
systematics at scission- 

Medium and heavy fissioning nuclei are deformed at their 
scission points and have less Coulomb energy than if they were i 
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FIGURE 3-1 Reduction of average fission-fragment kinetic 
energies by surface-plus-window dissipation, compared to 
experimental values. 




spherical- But simultaneously, their saddle shapes become less 

deformed and nucleus might accumulate prescission kinetic energy- 

The prescission kinetic energy must be smaller for viscous fluids 

as the amount of energy that is dissipated into heat is very 

> 

crucially dependent on the nautre and magnitude of viscosity 
prevalent in the process- In this way, estimation of nuclear 
viscosity can be made by systematics of total kinetic energy of 
symmetrically fissioning nuclei as a function of excitation energy 
ofthefissioningnucleus- 

The degree of dissipation of collective energy into excitation 
energy is also important in another context- Though the principal 
contribution to the kinetic energy comes from mutual coulomb 
repulsion, the collective energy of the vibrational states can also 
be transferred to the kinetic energy of the fission fragments. The 
translational energy of the descent of nucleus can also- contribute 
to the TKE - The single particle states also affect the final TKE- 
The contribution due to each can vary substantially depending on 
■the assumption regarding viscosity. Problem of separating from 
total kinetic energy of fragments all their component parts is very 


difficult. The 

experimental study for 

the 

effect 

of t ran s i t i on 

states of the 

fissioning nucleus on 

the 

total 

kinetic energy 


consists of the study of correlation between the mean total kinetic 
energy of the fragments and other characteristics of the fission 
process in which the properties of the transition state are 
reflected. 

The variation of the total kinetic energy as a function of 
excitation energy is available for a few fissioning systems- For 



example , iri the case of Pu , we have the data f or spontaneous 
fission and isomeric fission, and neutron induced fission of 
provide systematics upto 10 MeV of excitation energy- In isomer and 
resonance fission the excitation energy shows up entirely as an 

• I 

increase in kinetic energy- This led Swiatecki to 

the conclusion that such a fission process must be non—viscous • For 
neutron induced fission with excitation energies at or above 4-5 
HeV there is a sharp drop in kinetic energy by about 2-5 MeV 
followed by a gentle decrease- By simultaneous measurements of 
neutron energy with the assumption that average y energy remains 
the same, it was found that average internal excitation energy is 
increased by 5-4 MeV/, whereas kinetic energy by 1-15 MeV/- This type 
of fission process is highly damped- The sharp drop in kinetic 
energy at the excitation energy of 4-5 MeV/ can be attributed to the 
breaking of the nucleon pairs- The decrease in the average kinetic 
energy with excitation energy might be attributed either to the 
washing out of shell effect or to larger viscosity- Since breaking 
of nucleon pair is responsible for. transition from adiabatic to 
damped fission, adiabatic fission process should show odd-even 

2S5 

effect in mass yield and kinetic energy distributions- Wahl etal 
and Amiel and Feldstein^** measured the abundance of even Z products 
in thermal neutron fission of U and found it to be 20-505( larger 
than the relative average yield- Unik*^ found that fine structure 

' i . 

vanishes as we go to heavier nuclei *'***Cm to *'***Fm. David etal**”** 
have made extensive measurements on the variations of total kinetic 
energy as a function of excitation energy of the fissioning nuclei. 
They have employed the <a,a’f) reaction on several isotopes of 



of Uranium 


Uranium and Thorium to study this variation- In case 
nucleus, they find that the slope of the TkE' versus E curve is a 

X 

small negative number- The bar on the total kinetic energy shows 
the average value- For Th they find that TRF as a function of 

t 

excitation energy has a positive slope- David has also investigated 
the nature of the slope for mass separated fragments- They found 
that the slope of the TR£' versus curve is independent of the 
mass of the fission fragments- Hence, they conclude that the 
difference in the sign of the slope mentioned above is a property 
of the fissioning system- They suggest that this may have to do 
with the relative difference in height between inner and outer 
barriers. Back etal. measured the variation of total kinetic 

23S 

energy as a function of excitation energy for U in the 
excitation energy range of 5 to 50 MeV- Their low energy results 
corroborate with that of David’s and at energies above 22 MeV they 
find a small positive slope for the TKC versus curve- Choudhury 
etal-^"* have also measured the slope of the TR!F versus E^ curve in 
Pu- They studied the reaction <a,f) and not <a,a*,f> as that 
studied by David and Back- Moreover, they have measured the slope 
at energies much above the barrier- They also find a small positive 
slope for the curve- 

Although there is a plethora of experimental results available 
still it is not clear as to how the system becomes more viscous 
when it is given more excitation energy. Another open question is 
at what stage the system passes from a non-viscous fluid to a 
viscous fluid- If the total kinetic energy decreases with increase 
in the excitation energy, then how the system becomes more viscous 



as it is given more excitation energy ? If it increases with 
excitation energy then to what degree viscosity is playing a role 
in the transfer of energy from collective modes to single particle 
degrees of freedom ? And finally, why for some fissioning systems 

t 

total kinetic energy increases with excitation energy and for some 
it shows a decrease ? It is with the aim of trying to find an 
answer to some of these problems that the present study has been 
undertaken • In the present work we have measured the TRf as a 
function of excitation energy for for mass 5 ep«arated 

fragments- 

13* 21 EXPERIMENTAL DETAILS 

The experiment was performed at the Variable Energy Cyclotron, 
B«A-R-C«, Calcutta. The projectiles were 60 MeV a particles. The 

2 238 2 

target was 200 pg/cm U made by evaporating UF^ on a 100 /L/g/cm 
carbon backing- The experimental setup was mounted in the 60” 
scattering chamber in cave II of VECC- The chamber consisted of two 
movable arms each of which could be rotated by. 180*^ both from 
inside the cave as well as from the counting room- Schematic 
arrangement of the experimental geometry is shown in Fig. [3-23 • 

The telescope, consisting of a AE and an E solid state surface 
barrier detector, was used for measuring the energy of the 
inelast i cal ly scattered a particles and also for identifying 

theparticles- The AE detector was a fully depleted transmission 

* 

type silicon surface barrier detector while the E detector was a 3 
m-m- Silicon lithium drifted detector. Both the detectors were made 
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Fig. 3.2. Geometry of the experimental arrangement for the study of 
variation in total kinetic energy as a function of excitation 
energy of 238y 



at the detector laboratory of the VECC. The solid angle subtended 

by the telescope was 5 msr- The resolution of the E—AE particle 

telescope was about 1 MeV for the ine last i cal ly scattered alpha 

particles from the carbon backing- 

Two heavy ion surface barrier detectors F arid F were placed 

1 2 

as shown in the figure to detect the fission fragment pairs- They 

were placed at 180** to each other and each of them was 

approximately at 90^ to the beam axis- One of the detectors i-e- 

was placed very close to the target such that it caught fragments 

complimentary to those detected by F irrespective of the non 

2 

colinearity introduced due to the centre of mass motion and neutron 
evaporation from nascent fragments. The solid angle subtended by 
the fission detector F was 33 msr- The solid angle subtended by F 
fission detector was 14 msr- 

The angular distribution for inelastic scattering of an 

j 

energetic ion on a target material, due to incomplete momentum 
transfer, exhibits a <l-~Co59) behaviour- The detection of these 
ine las t i cal ly scattered alpha particles, if done at larger cross 
sections would help in increasing the count rate- Higher count 
rates would also result in a better true to chance ratio- Thus a 
proper positioning of the telescope is very important- For this a 
DWBA calculation was performed to calculate the cross section for 

230 

the ine las t i cal ly scattered alpha particles from U- It was found 
that the angular distribution is peaked at forward angles. However, 
a competing process, which contributes significantly to the 
deterioration of the true to chance ratio, is the inelastic 
scattering of alpha particles from the carbon backing which is also 



peaked at forward 

angles • 

Thus 

there 

is a trade off involved 

between higher 

count 

rates 

and 

contributions to 

chance 

coincidences- The 

optimum 

angle 

was fixe 

d at 30*" by actually 

taking 


the ratio of the two af oremerit ioned quantities during the actual 
2 

run- A 7 mg/ cm thick alumunium foil was mounted in front of the 
detector telescope to prevent the radiation damage of AE detector 
due to fission fragments reaching the detector- Magnets were also 
placed in front of the detector to suppress low energy electrons 
from the target- 

The electronic block diagram used for the experiment is shown 
in Fig-C3-3]- Fast coincidence was taken between the pulse from the 
fission detector which was further away and the AE detector while 
the geometry ensured the coincidence between E stop detector and 
nearer fission detector- The coincidence between the four pulses 
was further ensured software wise at the time of analysis- 

For the fast pulse, timing output from the preamplifier was 
fed to timing filter amplifier- The output from this was fed to a 
constant fraction discriminator in external delay mode- Walk 
adjustments were made to get a good timing resolution- The pulse 
from the preamplifier in the AE channel was processed in a similar 
fashion- An external delay was put in the AE circuit to bring the 
true TAC peak in the middle region of the TAC spectrum- The fission 
pulse started the time to amplitude coriverter and the AE pulse 
provided the stop pulse- 

Fig- [3-4] shows the TAC spectrum- The second peak is due to 
both the true and the random events whereas the other peaks are due 
to random events only- A random event is one in which an 
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ine las t i cal ly scattered alpha particle from carbon is detected in 
coincidence with fission resulting from the <:o(,f> process- As can 
be seen from the figure, the cross section for this kind of event 
is much larger than the events resulting from the <ct,a’f> process- 

The energy pulses from the two fission detectors, 
from the telescope E— AE, was fed to spectroscopy amplifiers- The 

i 

pulses from the amplifiers were fed to analog to digital converters 
(ADC) which were of the Uillkinson type, requiring a maximum analog 
signal of +8 volts- 

The TAC also provided a digital gate which was used as the 
master pulse which gated all the ADC’s and started the data 
acquisition cycle- Five pulses F^, F^, E, AE and TAC were acquired 
in list mode on a magnetic tape using a camac based data 

acquisition system- 

The experiment was continuously monitored on line by building 
(i>5 singles spectra namely F^, F^, AE, E and TAC Cii)3 dual 

spectra namely: 

(a)between F and F : to monitor the correlated fission 

12 

events- 

<b>between E and AE : to monitor the particle 

identification spectrum- 

(c>between AE and TAC s to monitor the true to chance ratio 

in the experiment. 

One conditional TAC spectrum was also built by cutting the 

elastic peak from the E spectrum to see the real true to chance 
ratio after eliminating the elastically scattered alpha particle. 

The beam current was kept at 5—10 t)A ; to avoid pile up in 



was 


fission as well as the AE detectors- The rate of true events 
around 8—10 per minute. 

[3.31 ANALYSIS PROCEDURE 

Mass spectrum from the correlated pulse height spectrum of the 

i 

two fission fragment detectors is obtained after calibrating the 
detectors and correcting the energies for : 

<a> energy loss suffered by the fragments due to finite target 
and backing thickness. 

<b) the centre of mass motion of the fissioning nucleus- 
<c) neutron evaporation from primary fragments- 

The analysis is involved as each of these corrections and the 
calibration depends on the mass of the fragment, something which we 
have set out to determine- We invoke the momentum and mass 
conservation to determine the mass of the fission fragments, and 
then apply these corrections iteratively till the mass of the 
fragments comes out to be the same in subsequent iterations- 
Appendix C2] gives the details of the analysis programme- It was 
found that after two iterations the value of the mass converges- 

E-3-3] <a) Energy Calibration! 

The output pulse height due to the energy lost by an energetic 
charged particle in a detector medium is a function of the particle 
energy, its mass and charge- For the lighter ions, the pulse 
produced in solid state surface barrier detectors is very nearly 



proport i onal to the energy of the particle but in the case of 
fission fragments and other heavy ions, particles having same 
energy but different masses do not produce the same pulse height 
and thus the linearity between the pulse height and the energy is 
lost- The pulse height produced in the solid state detector due to 
lighter ion is higher than due to a heavier ion of the same energy- 
The departure in linearity due to this dependence is known as the 
pulse h&ighi defect • This has to be taken into account in the 
energy calibration- 

It was found that once the mass is fixed the relationship 
between the pulse height and energy was unique- Further, once the 
energy is fixed, the relationship between pulse height and mass is 
linear- This was demonstrated by Schmitt etal^^ using Br and Iodine 
ions- Based on the above facts, there are a number of' calibration 
schemes, all empirical, but the one due to Schmitt etal has found 
the most widespread use- According to this procedure the 
calibration equation for a surface barrier detector reads as 

E = (a ^ a’M) X 4- b + b’M 


where X is the pulse height, H is the mass of the heavy ion, 
a,a’,b,b’ are the constants which have to be determined for each 
individual detector- Schmitt calibration provides a scheme by 
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which, measuring the fission fragment spectrum of U or Cf , 

one can calibrate their detector, if we know P , the pulse height 
of the most probable light fragments and P , pulse height of the 
most probable heavy fragments. The constants are given by the 
following relationsi 

a = a /<P - P„) 

O JL H 



a’= a’/CP - P > 

O L H 

b = b — a - P 
o L 

b’ = b» -a’ -P 
O L 

i 

where a , a’, b and b’ are constants given by Schmitt whicf' 
depend on whether one uses or for calibration. 

Though this scheme is most widely used, the values of the 
total kinetic energy obtained by this energy calibration scheme iT 
the double energy measurement and the one obtained due to the time 
of flight technique (i-e- the double velocity measurement), an 
known to be slightly at variance. To quote one result^^; th< 
average total kinetic energy release <TKE > for thermal neutroj 

235 

induced fission of U as obtained with surface barrier detector* 

^ . 

using this technique was reported to be < TKE > = 171-9 MeV-Froii 

double velocity measurement the figure turns out to be 168-3 MeVf 
Moreover, the calibration scheme of Schmitt was formulated in earl) 
sixties- Since then the detector technology has advanced much an< 
the material of Silicon may be different- Another drawback of thei? 
scheme is that they tested it for Bromine and Iodine ions only- 

Weissenberger etal^^ have recently reinvestigated thi 
behaviour of solid state surface barrier detectors with ORTEC I 
series heavy ion detectors- They have also checked for ions othei 
than iodine and bromine, and have arrived at the conclusion tha^ 
the phenomenological ansatz by Schimitt etal for the pulse heigh 
response of surface barrier detectors to heavy ions covering th» 
energy and mass range of fission fragments is perfectly wel 
suited. They conclude that within the experimental and evaluationa 
uncertainty of ±150 keV, there is no need to go beyond the linea 



parame tr i zat i on of the response E<X,M>- They have slightly modified 
the constants a^,a^’,b^ and b^’ for the pulse height response of 
Californium fission fragments- We have used these set' of constants 
in calibrating our detector: 

a = 24 - 30 
o 

a ’ = 0-0283 


b = 90-397 
o 

b ’ = 0-1150 
o 

For finding out P and P for the two detectors 

L H 

calibration was taken separately for the two detectors 
Guassians were fitted to calculate the position of P^ and 
peak positions for the light and heavy fragments. 

The Calibration of the telescope detectors was done by taking 
the energy spectrum of the inelatically scattered alpha particles 
on carbon nuclei(*^C). Since the energy level scheme of this is 
known, the difference in the peaks of the ine las t i cal ly scattered a 
particle spectrum is known, the calibration of the telescope could 

be done- 


a 

2=Zcf 

and 

then 


the 


[3-3] <b) Energy Loss Correction : 

The energy lost by a charged particle in a material medium is 
a function of its energy, mass and charge- Fission fragments with 
large A and Z and energy of the order of 1 MeV/nucleon lose a 
substantial amount of energy in the finite thickness of the target 
and the backing material- Apart from the energy loss, there is also 
energy straggling which is also dependent on the z of the material. 
This energy spread degrades the energy resolution which, in turn. 
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spoils the mass resolution- Thus in these experiments, where one 
measures mass distribution by double energy measurements in a 
typical back to back detector arrangement, the target and the 
backing should be thin and backing should be of low Z material- 
Initially, when the fragments are emitted they are fully ionized 
but as they transverse the medium they acquire electrons and lose 
their charge- Since the energy loss is a crucial function of the 
charge, in calculating it the term effective charge is used- 

There is a multitude of empirical relations for the stopping 
power calculations but unfortunately, there is a dearth of 
experimental results and most of the tables extrapolate into the 
regions uncovered by experiment. The rather extensive tables of 
NorthCliffe and Schilling present their data with unstated 

accuracy - Moreover , since their presentation, our knowledge of 
energy loss mechanism has improved. Hence, we have used the 

mechanism of Zeigler etal- , who have presented an extensive 

calculation for energy loss suffered by ions with different Z in 
different elements. For solid targets and ion velocities from 0-2 
to 2 MeV/u the stated error, at best, in their formulation is lOX- 

•K- 4f- 

The expression that they have used for Z^^ /Z^ is 

Z */Z * = 1 - exp<-A>C 1.034 - 0.1777 exp<-0.081 14>Z ] 

HI H ^ 

where 

A = B + 0.0387 sin<nB/2) 

B = 0.886 (E/25M^>‘^^/Z^*''®. 

CHere E is energy in keV, is mass in a-m.u. and is the 

charge. The expression for is 

(Stopping)"^ 



A 


for E/M < 999 keV/Za-m-u- 

Hhere S = A 

lov 2 

®high = > In [ 1 + A^E + A^/E ] 

Stopping = A //3^[ In ( ) _ ^ <lnE/ ] 

o 

The coefficients given for different elements. 

) E is the energy of the ion traversing the medium 
Z is the charge of the medium 

ft — v/c where v is the velocity of ion and c is the 
velocity of light. 

[3.33 < c>Centre Of liass Motion : 

The cx particle that hits the nucleus imparts some energy 

and momentum to it and gets inelasti cal ly scattered- Then, if the 
moving nucleus has energy above the fission barrier, it breaks into 
two fission fragments- The energy of the fission fragments that we 
measure is in the lab frame- We are interested in the. energy of the 
fragments in the centre of mass frame- To convert the lab energy to 
the centre of mass energy a small correction is required. Fig 3-4. 
illustrates how the centre of mass motion was taken into account- 
It is taken from Ref.[4.03 where such an iterative calculation for 
the mass of the fission fragments from double energy measurements 
has been described- If the compound nucleus was stationary the 
fragments would have been emitted at 180* with respect to each 
other, but since it is moving, they are emitted at an angle- The 
angle formed by the direction of their motion in the lab frame is 
known as the folding angle- 

Since we know the energy of the scattered a particle, thereby 




the energy imparted to the fissioning nuclei, we can calculate the 
velocity by formula 


V = 1 . 389 

GH 


CN GN 


Since we have measured the E and E the fission fragiiient 

1 2 

kinetic energies, we can calculate v and v using the above 

12 

relation between velocity, energy and mass- 

The velocity of the fragments , then, in the centre of mass 
frame is 


/ 2 2 
/ lab . Lab 
V + V 

1 2 


2v V CO50 > 
1 GN 1 


/lab 2 

/ V V 

1 GN 


2 V V cosO > 
1 CN 2 


It is sufficient to know the angle of emission of one of the 
fission fragments- The angle of emission of the other fragment can 
be calculated by simple vector algebra- Knowing the velocities one 


transfer them back to energy to get the centre of mass energy* 


[3-3] ( c> Neutron Correction : 

The major fraction of the energy released, not appearing as 
kinetic energy of the fragments, is dissipated by the emission of 
neutrons- The neutron, anisotropy studies have established that most 
of the neutrons are evaporated from fully accelerated fission 
fragments- The neutrons that are emitted from the fission fragments 
modify the energy and mass of the fragments- Since we are 
interested in the mass and energy distributions before neutron 
evaporation, we have to make corrections for the energy and mass as 
measured by our detectors for the energy taken away by the emitted 
neutrons- This correction is applied on the basis of previously 



measured neutron numbers* 

The number of neutrons that are emitted from the fission 
fragments is a function of excitation energy of the target nuclei 
and the kinetic energy of the fission fragments- At low energies 
i-e- energies near the fission barrier, number of neutrons is a non 
linear function of fragment mass- The average neutron yields are 
non integral implying thereby some variation in neutron yield from 
one fission event to another fission event- 

At low excitation energy of the compound nucleus, the most 
interesting aspect of neutron emission is the variation in neutron 
yield as a function of fragment mass- The first evidence for 
structure in the dependence of neutron yield on fragment mass was 
obtained by Fraser and Milton^^. There are unequal number of 
neutrons emitted from two complimentary fragments. Fig.C3.6<a>] 
(taken from Ref -[42] shows the variation of the neutron yield as a 
functioii of fragmerit mass for two different excitation energies for 
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the fission of U- This behaviour is exhibited by almost all 
fission fragments resulting from the fission of all actinides- This 

*43 

led Terrel to the conclusion that the neutron yields are closely 
related to def ormabi 1 ities of the nascent fragment- Closed shell or 
near closed shell nuclei F<refer spherical shapes and are unusually 
resistant to deformation. If the complimentary fragment of such a 
nuclei is not close to a closed shell, most of the deformation 
energy will be in the complimentary fragment. The number of 
neutrons depends on the amount of deformation energy stored at 
scission. Once the energy is increased, shell effects are 
diminished or get washed away and the closed shell nuclei are also 
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Fragment mass number 

Fig. 3.6(a). Variation of neutron yield as a function 
of fragment mass in (p + 2^®u) fission for 
proton energies of 11.5 MeV and 22 MeV. 



Fragment mass number 

Fig. 3.6(b). Variation of neutron energy with fragment 
mass. 
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able to share deformation energy, Fig- C3.6(a>3 shows that the saw 
tooth nature of neutron yield gets diminished as the excitation 
energy of the fissioning nuclei is increased- In our calculation, 
we require i^<A, where u is the number of neutrons, A is the 

fragment mass and is the excitation energy- Unf ortunately , there 
is a dearth of experimenal results and there is no emperical 
relation to give the number of neutrons as a function of fragment 
mass and excitation energy of the fissioning nuclei- We have 
therefore, extrapolated or in trapo lated linearly the results of 
obtained by Bishop etal depending on the excitation energy of the 
Uranium nuclei by fitting least square fits- 

The average energy taken by the neutrons in the lab frame is 
about 2 MeV- Since most of the neutrons are emitted from moving 
fission fragments, the energy given to them by virtue of fragment 
motion is 2/3 MeV- Fig-[3-6(b)] shows the average neutron energy as 
a function of fragment mass* The data is again from Bishop etal 
The least square fit to this data has a very small slope* Though 
within the experimental resolution and the evaluational 
uncertainties we could have assumed that the energy is constant, 
yet we have taken a small slope* The energy is given by the 
relation 

C 50 ) X 4*41 17647* 10~* + 1*1268382 

where is the mass of the fragment. 

13.41 RESULTS AND DISCUSSION 


Fig.C3-7] shows the variation of average total kinetic energy 




Excitation energy , MeV 


Fig. 3.7. The variation of total kinetic energy with excitation energy. 
The o show the experimental data , the X shows the 
data corrected for neutron from saddle to scission , the 
A shows the data further corrected for scission neutrons. 
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(TRF>as a function of the excitation energy of the fissioning 
nucleus- We observe that the total kinetic energy in the excitation 
energy region of 8-5 to 12-5 shows a decrease with excitation 
energy- Here, we have neglected the point at the excitation energy 
of 6-5 HeV because the error due to chance coincidence is very 
large due to the inelastic peak from the carbon in the backing 
interfering at this value of the excitation energy- This implies 
that as the excitation energy is increased, the viscosity of the 
system increases- This has been observed by Back etal^^ and also by 

2lP 

David etal- - The total kinetic energy then remains constant till 
22-5 MeV and then it shows an increase with the excitation energy- 
At high excitation energies Back etal- also observe a positive 
slope- The positive slope has also been measured by Choudhury 
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etal , though in a different system and at slightly higher 
excitation energies of the fissioning system- A possible 
explanation for this change of slope as more excitation energy is 
given to the system could that at lower excitation energies the 
shell effects are prevalent and giving a small amount of excitation 
energy to the system may amount to breaking of paired nucTeons- At 
22 Mev the pairing correlation is destroyed and the shell effects 
get diluted- They are not yet completely wiped out as is evident 
from the fact that the mass distribution still shows two humps- The 
peak to valley ratio averaged over all the values of excitation 
energy is 1 is to 1.5- The dilution of shell effects may weaken the 
coupling of the internal degrees of freedom to the external 
degrees. The excitation energy of the fissioning system then may 
manifest itself in the relative kinetic energy of the fragments. 
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Table -3-1 

Tab 1 e - 3 - 1 
f urict i ori 
Fig. [3-8] 
var iat i on 
f ragment s 


The variation of the total kinetic energy as a function 
of excitation energy for mass selected fragments- 

shows the variation of the total kinetic energy as a 

of excitation energy for mass selected fragments- 

shows the data of Table [3-1]- The figure shows the 

of TRF with excitation energy for mass selected 

For this kind of analysis the statistics of our data is 


very poor and consequently the error bars are very large 



























































































Nevertheless, within the experimental errors one can say that for 
mass selected fragments near the doubly closed shell the total 
kinetic energy does not show any variation with excitation energy® 
For the highly asymmetric masses, the total kinetic energy 
increases with excitation energy of the fissioning nuclei* This 
implies that the coupling between the internal and the external 
degrees of freedom is stronger in case of spherical miclei as 
compared to deformed nuclei- 

Chance coincidence correction 

2^***33 

One drawback of the earlier measurements is the 
corrections applied to eliminate the effect due to chance 
coincidences. These chance coincidences, as mentioned earlier, 
arise due to alpha scattered from the carbon in the backing and 
fission due to <a,f> process- In earlier measurements, correction 
for the chance coincidence is made by subtracting the events under 
one of the random peaks from the events under the true + random 
peak after processing both the events as mentioned in Sec- [3-33 • 
This is improper because events under the random peak and random 
events in the true + random peak are due to the (a,f> process where 
the centre of mass and the neutron corrections apply differently 
than in the <a,a’'f) process- The random events under the true + 
random peaks have therfore to be processed in the same manner as 
those events under the random peak- This is not possible to do- In 
the earlier measuremeri ts , the correction for this effect has been 
done assuming that the data is the result of the direct reaction. 
This subtract ion results , therefore, in biasing the data- Even if 
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the true to chance ratio is one is to ten, it introduces an error 
of 10 % in their results. We have therefore not done this 
subtraction as this results in undue biasing of the data- 

Thus we suggest that this expter imen t should be done with a 
self supporting foil of the target nucleus- This has not been done 
earlier because the target would have to be much thicker- Being a 
high Z material, the fragments will loose a substantial portion of 
their energy in the foil. This again will increase the uncertainity 
in the results. However, if one uses a gas detector to measure the 
energy and also the angle of emission of the fragments, the error 
will then be only due to the energy straggling which is of much 
smaller magnitude- The measurement of angle will also result in 
substantial improvement of the mass resolution since the averaging 
over the angle, which is done to correct for the energy loss 
correction and centre of mass motion in the experiments so 
far, spoils the mass resolution- If this experiment is done with 
mass resolution of the order of one mass unit, it may be possible 
to observe the effect of pair breaking on the total kinetic energy 
- an effect which has not been seen earlier due to poor mass 
resolution- 

Neutron Cor recti on: 

We have made corrections for the number of neutrons emitted by 
the fully accelerated fragments- The number of neutrons emitted by 
the fragments at low excitation energy of the target is a function 
of their mass as well as the excitation energy of the fissioning 
nuclei- These neutrons can be classified into three categories 



depending on the (tie chan ism of their- emissions 
(i> neutrons emitted from the compound neucleus 

(ii) neutrons emitted when the fissioning nucleus descends from 

saddle to scission 

< i i i )neutror!5 emitted from fully accelerated fragments 

The study of prompt neutrons and their angular distributions 
in early sixties^^ had revealed that a small fraction of 

neutrons do not come’ from the fragments, but are emitted by the 
fissioning nucleus by the first two mechanisms mentioned above- The 
neutrons emitted before scission are nearly isotropic in the c-m 
system of the compound nucleus, whereas neutrons emitted by fully 
accelerated fragments are strongly correlated to the direction of 
fragment mot ion - 

The measured fragment energies have to be corrected for 
neutron emission because the neutrons are emitted with some kinetic 
energy and in addition to this have some energy due to emission 
from fragment motion- 

The results are corrected using the experimental results of 
Bishop etal ■ They have measured the number of neutrons emitted in 
the proton induced fission of Uranium for two different proton 
energies- We have interpolated or extrapolated for the number of 
neutrons basing our results for these two energies following a 
procedure very similar to the one used by David - Unfortunately 
the data of Bishop etal does not distinguish between the neutrons 
emitted before scission and the ones evaporated from fully 
accelerated fragments- Neutrons emitted from fully accelerated 
fragments are to be taken into account for the correction of the 
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single fragment kinetic energies- This experimental data does not 
exist- So far all the existing results on TKE measurement s suffer 
from this uncertain i ty of making a correction to single fragment of 
the type ( 1+ v>/M > (where u is the number of neutrons, li is the 
mass and is the energy of the fragment which is overestimated 
because it includes contributions from mechanisms 1 and 2 which do 
not modify the single fragment kinetic energy- 

To correct for the number of neutrons emitted by the 
fissioning nucleus from compound nucleus to the saddle point stage, 
we have made a calculation of number of neutrons emitted from 
compound nucleus using F and F values from Wandenbos ch"^^. The 

n f 

computer code to calculate this is given in Appendix [3]- If we 
incorporate this calculation in our experimental curve showing 
TKE vs E , the curve which was having a positive slope becomes 
nearly a constant- The value of the TKE as a function of excitation 
energy, after correcting for this effect, is shown in Fig-C3-7] 
with a X- Our calculations do not include neutrons emitted from the 
saddle to scission- 

There is a lot of data in heavy ion p»hysics on prescission 
neut rons'^^'*^- Gavron etal have found that theory always 
underestimates the value and most of the excitation energy goes 
into the emission of prescission neutrons. There is, however, one 

' 50 ' 

experimental result in a induced fission by Scobel etal, which has 
direct bearing on our experimental results- They have studied 
prescission neutrons from fission fragments from a induced fission 
of Uranium. They find that nearly half of the excitation energy 
appears in the form of enhancement in prescission neutrons- We have 
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^gaiii corrected our results based on the results of Scobel which 
also gives the average energy taken by the scission neutrons as a 
function of excitation energy. This we show in the Fig- [3-7] with 
A* However, they have not measured the average energy and the 
multiplicity in low excitation energy region- Thus to quote the 
value of slope for Tk E / E one requires the measurement of the 

X 

prescission neutrons in this excitation energy range- For depending 
on the multiplicity of these neutrons the curve can have a 
negative, positive or zero slope- 
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CHAPTER IV 


SUMMARY AND CONCLUSIONS 


The studies in this thesis were aimed at providing further 
insight into the fission dynamics during the descent of the 
fissioning nucleus from the saddle to the scission point- We have 
investigated the prompt gamma ray emission in spontaneous fission 
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of Cf and the total kinetic energy distribution in the fission 
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of U following inelastic alpha scattering. 

The results obtained can be summarized as follows- 
(i) The gamma ray multiplicity distribution measurements indicate 
that the spin distribution of the fission fragment spin has a width 
(variance) of <7 = 5-7 ± 0-5 is deduced from the data for the charge 
ratio 1-3 (heavy fragment charge=56 and light fragment charge==42>» 
This is a new piece of information since all the earlier studies 
have yielded only the mean value of the fragment spin- 

The value of the width obtained is in general agjreement with 
the hypothesis (Moretto etal«^> that the different angular momentuin 
bearing modes (namely (Dtwisting (ii> bending (iii) tilting and 
(iv> wriggling) at scission are in equilibrium- This width has been 
compared with the eKpression given by Moretto assuining t|ie scission 
temperature to be around 0-5 to l-O He\/- 



Gamma ray angular distribution shows that the anisotropy with 
respect to the fragment direction has a strong quadrupole signature 
in the range of 300<E^<700 keV the maximum being l-26±0-01 around 
500 keV- Beyond 700 keV the anisotropy decreases but remains more 
than one- Below 300 keV the anisotropy is less than one signifying 
a predominantly dipole emission- 

The analysis of the gamma ray data has brought out a number of 
shortcomings of the existing statistical model calculations of 
gamma ray emission- This may have important consequences on the 
de—exci tation mechanism of nuclei having a finite spin and 
excitation energy- 

The anisotropy of the prompt gamma rays with respect to the 
fragment direction for gamma energies greater than 700 keV being 
more than one is not understood- The statistical model, with the 
avai lable ^ expressions for the decay widths of gamma emission, 
predicts that the emission is predominantly dipole in nature at 
energies beyond 700 keV- This should make the anisotropy less than 
one- The experimental data suggests about an equal admixture of 
dipole and quadrupole transitions with the probability of 
quadrupole being slightly more- The origin of these quarupole and 
their nature is presently not known - 

Deformed nuclei and spherical nuclei having large excitation 
energy and spin differ in their deexcitation- The present 
statistical model simulation does not include any difference in the 
nature of de— excitation between the spherical and deformed nuclei- 
Also the fragments carry angular momentum both as aligned single 
particle spins and collective rotations <in case of deformed 



f ragmen ts > - The answer to how a given angular momentum is 
partitioned between the two is also not answered by the theory- 

Investigations on the kinetic energy distribution of the 
fission fragments have brought out the following resultss 

The total kinetic energy in the excitation energy range of S-5 
HeV to 12<.5 MeV shows a decrease with excitation energy- In the 
excitation energy range of 12-5 MeV to 22-5 MeV it remains 
constant- From 22-5 lieV to 28-5 MeV it shows an increase with 
e X c i tat i on energy - 

The excitation energ9 dependence of the mass distribution 
leads to an apparent dependence of the fragment kinetic energy on 
the excitation energy- To minimize this effect, we carried out the 
variation of total kinetic energy as a function of excitation 
energy for mass separated fragments- We found that for nuclei near 
the doubly closed shell, the total kinetic energy does not show as 
much variation with excitation energy as it does for more 
asymmetric masses- While our statistics for a complete detailed 
analysis of this kind was poor, within experimental errors we can 
say that the descent of the fissioning nucleus from the saddle to 

scission for near spherical fragments is more viscous than for more 

asymmetric fragments- 

Further analysis has brought out the reliability of earlier 

i , 

results*"* in question. This relates to the correction to single 

fragment kinetic energy due to neutron emission. The multi chance 

fission events do not modify the single fragment kinetic energy. 
However, these have been included in the correction for the single 
fragment kinetic energy in the previous results. 
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We have, unlike ihe previous measur ement-s , corrected our curve 
of Tlk£ versus excitation energy partly by excluding the neutrons 
emitted from compound nucleus to the saddle point- The correction 
was done on the basis of statistical decay widths for fission and 
neutrons- Once the neutrons emitted during this stage of the 
fission process are excluded from the neutron correction in the 
single fragment kinetic energy, the nature of the curve beyond 22-5 
MeV changes considerably- 

Another new dimension has been added by the presence of 
enhanced prescission neutron multiplicity observed in heavy ion; 
reactions, an area which has received very little attention in the 
past but plays a very crucial role in the analysis of this kind of 
measurements- These neutrons like those emitted in multichance 
events do not modify the single fragment kinetic energy- 

From the measurements of Scobel etal-**, we have corrected our 
total kinetic energy variation with excitation energy for errors 
resulting due to inclusion of these neutrons in the correction of 
the single fragment kinetic energy. This makes the curve of fl<iC/E^ 
almost flat- A reliable information on the viscosity prevalent 
during the descent of the fissioning nucleus from the saddle to the 
scission point can be deduced from these results, if the 
experimental information on the dependence of scission neutron 
multiplicity on the excitation energy of the fissioning nucleus and 
the fragment mass is known- 
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APPENDIX 1 


C PROGRAM TO CALCULATE THE SPIN DISTRIBUTION 

PROGRAM WIDTH 

DIMENSION SNEU<1000),SE1S<1000),SE2S<100>,SE2CC100>,RAND<100G> 
DIMENSION NITE1C15>,NITN<15>,TACC100),NSPEC<50>,NN(10>,NU<10) 
DIMENSION AZ<50> 

CHARACTER* 10 NAME 

COMMON/AA/HCROS, PI ,SQRPI ,PISQR, DELTA 

C0MM0N/BB/A,Z,EINI2,JINI2,ABAR,RINER,BN<3,50>,ITZ,INEU 
EXTERNAL FNEU, FEIS, FE2S, FE2C 

PRINT* <”*GIVE THE NAME OF THE OUTPUT FILE:”)’ 

READ ’<A10)*,NAME 

100 F0RMAT<A10> 

101 FORMAT <2X, ’MASS OF THE NUCLEUS=’ ,F5.0,2X, ’CHAR6E=’ ,F5.0,2X, 

1 ’EX. ENERGY=’,F<t.0,2X,’SPIN=’,I4,2X,’DELTA=’,F5.1> 

102 F0RMAT(5X,<;4C14<1H->),*X’>) 

103 F0RMAT<4X, ’INITIAL EX- ENERGY=’ , F8.3) 

104 F0RMAT<10X, ’SPIN VALUE IS=*,I3) 

105 F0RMAT(4X,’ THE SHELL CORRECTION PARAMETER=’ , F6.3) 

106 F0RMAT<4X, ’FIRST AND SECOND MOMENT OF THE MULTIPLICITY DISTRIBUT 
IION OF GAMMA RAYS ’> 

107 FORMAT <8X, ’FROM FISSION FRAGMENTS AT DIFFERENT EXCITATION ENERGY 
1 AND SPIN’) 

108 F0RMAT<5X, ’EX ENERGY’ ,5X , ’NEUTRON’ ,5X , ’STA E1’,5X,’STA E2’,5X, 

1 ’COL E2*) 

109 FORMAT <2X, ’NO. OF GAMMA RAYS=’ , 12, 2X ,’ ENERGY OF THE El STATISTICAL 

1= ’,F5.3) 

110 F0RMAT<2X, ’NO-OF NEUTRON=’ , 13, 3X ,’ ENERGY OF THE NEUTRON=’ ,F5.3 
1 ) 

0PEN<UNIT=4 , F I LE=NAME > 

0PEN<UNIT=33,FILE=’BIND:DAT’ ,ACCESS=’READ’) 

DO FOR 1=1 ,27 

READ<33,*)AZ(I),BN<1,I),BN<2,I),BN<3,I) 

ENDDO 

CL0SE<UNIT=33) 

WRITE<4,106) 

WRITE<4,107) 

WRITE<4,*) 

WRITE<4,*> 

IN=1 

PI=4.*ATAN<1.) 

PISQR=PI**2 

SQRPI=SQRT<PI) 

RMASS=1.04 

ITER=0 

PRINT* <”*GIVE A RANDOM NUMBER AN INTEGER :”)’ 

READ*, IX 

ITER=0 

Z=62.0 
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E INI 1=20-0 

EINI=EINI1 

JINI1=15 

JINI=JINI1 

IAZ=1 

A=AZ<IAZ> 

AINI=A 

DELTA=4.0 

WRITE<4, 101 )A,Z,EINI1,JINI1, DELTA 

HCR0S=6.5819 

IAZ*IAZ+1 

ITZ-1 

GO TO 18 

10 IF<Z.LE. 50-0)00 TO 60 
ITER=0 

Z=Z-1 -0 
A=AZ<IAZ> 

AINI=A 

IAZ=IAZ+1 

WRITE<4,*> 

WRITE<4,*> 

WRITE<:4,101>A,Z,EINI1,JINI1, DELTA 

ITZ=ITZ+1 

JINI=JINI1 

EINI=EINH 

C WRITEC4,108) 

DO 11 IS=1,10 

11 NN<IS>=0 
GO TO 18 

12 ITER=0 
EINI=EINI-5- 
WRITE<4, 103)EINI 
JINI=JINI1 

IF<EINI-LE-5-)G0 TO 10 
GO TO 18 

13 ITER=0 
JINI=JINI-5 
IF<JINI-LE-0)G0 TO 12 
WRITEC4,104)JINI 

C PRINT 164,JINI 

GO TO 18 

14 ITER=0 
MSUM=0 
NSUH=0 

DO 15 L=l,10 
NSUM=NSUM+NU< L ) * < L-1 ) 

15 MSUM=MSUM+NN(L>*(L-1) 
RSS=FLOAT<MSUM) 

SSS=FLOAT<NSUM) 

RSUM=RSS/100. 

SSUM=SSS/100- 

RTT=0-0 

STT=0-0 
DO 16 L=l,10 
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STT=STT+NUCL>*<SSUM-L+1 >**2 
ie> RTT=RTT+NN<L>*CRSUM-L+1)**2 

XTT=SQRT<RTT/100. ) 

YTT=SQRT<STT/ 100 . > 

WRITE<4,*> RSUM , X TT , SSUM , YTT 

DO 17 LL=1,10 

NU(LL>=0 

17 NN<LL>=0 

C PRINT ♦.DELTA 

DELTA=DELTA-3.0 
IF<DELTA.LT. .0>THEN 
DELTA=4.0 
GO TO 13 

C GO TO 60 

ENDIF 

18 ITER1=0 
A=A1NI 
INEU=1 

RADIUS=1 .2^A^^<1 . /3. > 

R I N ER= . 4^ A^RM ASS^ ( RAD I US*^2 )/hcros^^2 
HINDl=<10.^^<-3>) 

HIND2=<10.^^<-5>> 

ABAR=0.14^A 

GB AR=6 . ♦ AB AR / P I SQR 

ITER=ITER+1 

IF<ITER.GT.100>G0 TO 14 

EINI2=EINI 

JINI2=JINI 

NNEU=0 

NE1S=0 

NE2S=0 

NE2C=0 

19 RADIUS=1.2^A^*<l./3.) 

RINER= . 4^A^RMASS^ ( RAD I US^^2 > / h c r-o s^^2 
ABAR=0.14^A 
GB AR=6 . ♦ABAR/ P I SQR 
ER0T=JINI2^^2/<2.»RIiMER> 

IF<<E1NI2-ER0T-.5).LE.0.0>THEN 

NN<NE1S+1)=NN<NE1S+1)+1 

NIJ<NNEU+1>=NU<NNEU+1) + 1 

NE2C2=JINI2/2 

GO TO 18 

ENDIF 

T=SQRT<EINI2/ABAR) 

CE1S=5. 1^<A^^<2./3.>>*<GBAR*^2>^HIND1*<10-**C-10>> 

GAUN=GINTEG< FNEU , 0- , 2 • ♦T>+GINTEG<FNEU , 2 - ♦T , 4 - ♦T)+GINTEG< FNEU 
1 ,4.^T,20.^T) 

GAMN=GAUN/WEJ<EINI2,JINI2> 

6AUE1S=GINTEG<FE1S,0.,2.^T>+6INTEG<FE1S,2.^T,4-^T)+GINTEG<:FE1S 
1 ,4.^T,6.^T>+GINTEG(FE1S,6.^T,3.^T>+GINTEG<FE1S,8.^T,20.^T> 

GAME1S=GAUE1S/WEJ<EINI2,JINI2> 

GTOT AL=GAMN+GAME 1 S 

PWNEU=GAMN/GTOTAL 

PWE1S=GAME1S/GT0TAL 
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ITER1=ITER1+1 

ZZ=RNDM<:iX,IN> 

IFCZZ.GT.O.O.AND.ZZ-LE.PWNEU)GO TO 20 
IF<ZZ.GT.PWNEU.AND.ZZ.LE. <:PWE1S+PWNEU>>G0 TO 30 

20 CONTINUE 
SNEU(1>=0.0 

DO 21 1=2,1000 

EPS=I*0.01 

GAM=FNEU<EPS) 

21 SNEU(I>=SNEU<I-1>+GAM 
DO 22 1=2,1000 

22 SNEU<I)=SNEU<I>/SNEU<1000> 

ZZ=RNDM<IX,IN) 

ZZ2=ZZ 

DO 23 1=2,1000 
IFCSNEU<:i>.GE.2Z>G0 TO 24 

23 CONTINUE 
IF<I.GE-1000)THEN 
EN=EINI2 
NNEU=NNEU+1 
A=A-1.0 

GO TO 18 
END IF 

24 EN=I*0.01 
NNEU=NNEU+1 

EINI2=EINI2-EN-BN<INEU,ITZ) 

INEU=INEU+1 

A=A-1.0 

C WRITE<4, 110)NNEU,EN 

GO TO 19 

30 CONTINUE 
SE1S<1)=0.0 

DO 31 1=2,1000 

EPS=I*0.01 

GAM=FE1S<EPS) 

31 SE1SCI>=SE1S<I-1)+GAM 
DO 32 1=2,1000 

32 SE1S<I>=SE1S(I>/SE1S<1000) 

ZZ=RNDM<IX,IN> 

2Z3=ZZ 

DO 33 1=1,1000 
IF<SE1S(I>.GE.ZZ>G0 TO 34 

33 CONTINUE 

34 EiS=I*0.01 
NE1S=NE1S+1 
EINI2=EINI2-E1S 

C WRITEC4,109)NE1S,E1S 

GO TO 19 

40 CONTINUE 
SE2S<1)=0.0 
DO 41 1=2,100 
EPS=I*0-1 

GAM=(EPS**5)*WEJ<EINI2-BN<INEU,ITZ>-EPS,JINI2> 

41 SE2S<I>=SE2S(I-1>+GAM 
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42 

43 
45 


50 

51 

52 

53 

54 


60 


DO 42 1=2,100 
SE2S<I>=SE2S<I)/SE2S<100> 

ZZ=RNDM<IX,IN) 

ZZ4=ZZ 

DO 43 1=1,100 

IF<SE2S<I> .GE.ZZ>GO TO 45 

CONTINUE 

E2S=I*0.1 

IF<I.EQ.101>THEN 

E2S=0 . 1 

NE2S=NE2S+1 

WRITE<4,*) ’ENERGY OF THE E2 STATISTICAL ’, ITERl , E2S 
GO TO 18 
END IF 

NE2S=NE2S+1 

WRITE<;4,*) ’ENERGY OF THE E2 STATISTICAL ’, ITERl , E2S 

EINI2=EINI2-E2S 

GO TO 19 

CONTINUE 

SE2C<:i>=0.0 

DO 51 1=2,100 

EPS=I*0.1 

GAM=EPS**5 

SE2C< I )=SE2C< I-l >+GAM 
DO 52 1=2,100 
SE2C<I)=SE2C<I>/SE2C<100) 

ZZ=RNDMCIX,IN) 

DO 53 1=1,100 

IF<SE2C(I> .GE.ZZ)GO TO 54 

CONTINUE 

E2C=I*0. 1 

IF<I.EQ.101>THEN 

E2C=.l 

NE2C=NE2C+1 

WRITE<4,*) ’ENERGY OF THE E2 COLLECTIVE ’, ITERl , E2C 
GO TO 18 
END IF 

NE2C=NE2C+1 

WRITE<4,*) ’ENERGY OF THE E2 COLLECTIVE ’, ITERl , E2C 

EINI2=EINI2-E2C 

JINI2=JINI2-2 

GO TO 19 

CONTINUE 

STOP 

END 

FUNCTION GINTE6<FF,A,B) 

DOUBLE PRECISION XX, WX 
DIMENSION XX<40>,WX<40> 

DATA XX/O- 0950 125098D0,-0. 0950 125098D0, 0.2S16035508DO, 

1 -0.2816035508D0, 0.458016777700,-0- 4580 167777D0, 

2 . 0.6 17876244400,-0. 61 78762444D0, 0.7554044084D0, 

3 -0.7554044084D0, 0.865631202400,-0-865631202400, 

4 0. 944575023 IDO, -0.9445750231 DO, 0.9894009350D0, 

5 -0-9894009350DO, 24*0- DO/ 



<J1 W N) 


2 


C 


C 


DATA WX/0. 1894506 105D0, 0 - 1894506 105DO , 
1 0-1 826034 150D0, O- 1691565194D0, 

0 - 1 495959888D0 , 0-1 495959888D0 , 
O. 12462897 12D0, 0-09515851 17D0, 
0.0622535239D0, O - 0622535239D0 , 
0 . 027 1 524594D0 , 24*0 - DO/ 

T=0.0 

DO 2 1=1,16 

XI=(B-A>*XXCI)*0.5+<B+A>*0.5 


0- 1 826034 150DO, 
0.1691565194D0, 
0-1 2462897 12DO, 
0.09515351 17DO, 
0.0271524594DO, 


T=T+FF<XI)*WX<I> 

CONTINUE 

GINTEG=T*<B-A)*0.5 

RETURN 

END 

FUNCTION WEJCEXjJ) 

COMMON/AA/HCROS , PI , SQRPI , P ISQR , DELTA 

COMMON/BB/A , Z , E INI2 , J INI2, ABAR , RINER , BN<3 , 50 > , ITZ , INEU 
ER=J**2/ <2 - *RINER) 

EF=EX-ER 

WEJ=0.0 

IFCEF.LE.O.O)RETURN 

FABAR=ABAR*< 1-<DELTA/EF>*< l-EXP<-EF*0-054> ) ) 
WEJ=<SQRPI*EXP<2.*<<FABAR*EF>**0.5)>>/(12.*<FABAR**.25>* 

1 CEF**1.25)) 

RETURN 

END 

FUNCTION FNEU<X) 

DIMENSION YY(7) 

COMMON/BB/A, Z,EINI2,JINI2, ABAR, RINER, BN<3,50>, ITZ, INEU 
COMMON/AA/HCROS , PI , SQRPI , P ISQR , DELTA 
DATA YY/31. 05, -25-91, 342-4, 21 .89,0.223,0-673,617.4/ 
Y=EINI2-BN<INEU, ITZ>-X 

SIGN=YYC3>*<:<A>**<l/3>>+YY<4>*<<A)**<2/3>>+(YY<l>*<;<A)**C-l/ 
1 +YY<2>>*X+CYY<5>*<<A>**<4/3)>+YY<6>*<<A)**<2/3>>+YY<7)>/X 

FNEU=<1./<<PI*197.)**2>>*940-*X*WEJ<Y,JINI2>*SIGN 
FNEU=X*WEJ(Y,JINI2> 

RETURN 

END 

FUNCTION FE1S<X> 

COMMON/BB/A, Z,EINI2,JINI2, ABAR, RINER, BNC3, 50), ITZ, INEU 
COMMON/AA/HCROS , P I , SQRP I , P I SQR , DE L TA 
ERR=34 - *A** < - 1 /6 ) 

TA=6. 

SIGS=60.*<CA-Z)*Z/A>*<2./(PI*TA))*C<<TA*X>**2)/<CTA*X>**2+< 

1 <X)**2-CERR>**2)**2>) 

FElS=<;i./<PI*197.>**2>*SIGS*CX**2)*WEJCEINI2-X,JINI2> 

FE1S=((X>**4>*WEJ<EINI2-X,JINI2> 

RETURN 

END 

FUNCTION FE2S<X> ' 

COMMON/BB/A ,Z,EINI2,JINI2, ABAR , R I NER , BN < 3 , 50 ) , I TZ , INEU 
FE2S=(<X>**6)*WEJ<EINI2-X,JINI2> 

RETURN 

END 
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FUNCTION FE2C(X> 

C0MM0N/BB/A,Z,EINI2, JINI2,ABAR,RINER,BN(3,50) ,ITZ, INEU 

FE2C=<X>**5 

RETURN 

END 

FUNCTION RNDM<IX,IN> 

INTEGER*4 JM,MM,IX,IC 
C PRINT *,IX 

IF<IN> 30,40,30 
30 IN=0 

JM=2**10+3 
MM=2**< 20 > 

FM=MM 

40 IC=IX*JM/MM 

IX=IX*JM-IC*MM 

FX=IX 

XD=FX/FM 

RNDM=ABS<XD> 

C PRINT*,RNDM 

RETURN 
END 
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PROGRAM TAPETEST 
INTEGER*2 IBUF1<1024> 

INTEGER*4 IBUF2<80> 

DIMENSION NTKE<2, 20, 1550>,NM1<2, 20, 1200), IE1C2, 20,1200) 
DIMENSION RMEAN3<2,20),IE2C2,20,1200),IMEA<2,100,20) 

DIMENSION IMRE<2,20,100,100),IME<2,20,100,100) ,ITAC<:750) 
DIMENSION IF1<20,200),IF2(20,200),IF3(20,250),IF4<20,200) 
DIMENSION IFIS1<:2048),IFIS2<2048),IPAR<2048) 

COMMON/AA/EP 

COMMON/BB/SU<20) 

COMMON/ CC / SC < 20 ) 

COMMON /DD/ EX 

COMMON/ EE /THETA 1 ,THETA2 

COMMON/ EF /THETA 

C0MM0N/NEUT/ECM1,ECM2,E1N,E2N,RNU1 ,RNU2,AM1 ,AM2 
INTEGER*4 lERRNO 
CHARACTER*4 C0MMAND<6) , ANS 
CHARACTER* 1 RAND 

DATA COMMAND/ ’BKSP’ , ’RWND’ , ’READ’ , ’EXIT’ , ’SKIP’ , ’WRIT’/ 

IC0UNT=0 

JK=0 

OPEN <UNI T=5 , F I L E= ’ POT : SYMB ’ , ACCESS= ’ READ ’ , STATUS= ’ OLD ’ ) 
OPEN<UNIT=G , F I LE= ’ MASS : OUT ’ , ACCESS= ’WRITE’, STATUS= ’ UNKNOWN ’ ) 
OPEN < UN I T=7 , F I L E= ’ MAS9 : OUT ’ , ACCESS= ’WRITE’, STATUS= ’ UNKNOWN ’ ) 
OPEN< UN I T=8 , F I LE= ’ TKE 1 : OUT ’ , ACCESS= ’WRITE’, STATUS= ’ UNKNOWN ’ ) 
0PEN<UNIT=9,FILE=’PARTI :OUT’ , ACCESS= ’ WRITE ’, ST ATUS=’ UNKNOWN’ ) 
0PENCUNIT=10,FILE=’TKEM:0UT’ , ACCESS=’WRITE’ , STATUS= ’ UNKNOWN ’ > | 
MTUR=40B 

1000 FORMAT</, ’«GIVE FUNCTION BKSP/RWND/READ/EXIT/SKIP/WRIT ??’) 

1001 FORMATC’ Header Words Read: ’,16) 

1002 FORMAT<’ No. of Blocks Read:’, 15) 

1003 F0RMAT<120(1H-)) 

1004 F0RMAT<25X,’ VALUE OF TOTAL ENERGY=’,I2) 

1005 F0RMAT(20I4) 

1006 F0RMAT<30X,’ ENERGY=’,I5) 

1007 F0RMAT<2X , 15 ,5X , ’ SU.M=’ , F 10 .3 , 5X , ’MEAN=’ , F lO -3 ,5X , ’ SIGMA= ’ , F 10 

1008 FORMAT <35X, ’MASS =’,I5) 

1009 F0RMAT(7X,I5,5X, ’SUM=’ ,F10.3,5X, ’MEAN=’ , F 10.3,5X , ’SIGMA=’ , FIO.; 
READ<5,*)<SU<I),I=1,20) 

READC5,*)<SC<I),I=1,20) 

CLOSE (UNI T=5> i 

FM1=1174. iFM2=1204. ; 

C FM2=1150. iFMl=1170. ; 

AC1=1146.0 * 

BC1=1554.0 

AC2=1080.0 
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BC2=1496.0 

EP=60.0 

PI=4.0*ATAN<1.0> 

THETA 1=90-0 

THETA2=90.0 

EH 1=78- 42 

EH2=78.42 

EL1=102.64 

EL2=102.64 

TA1=<EL1-EH1>/<BC1-AC1> 

TA2= < E L 2-EH2 > / C BC2- AC2 ) 

TB1=EL1-TA1*BC1 

TB2=EL2-TA2*BC2 

AN0T=24.3 

AN0TP=0.0283 

BN0T=90.397 

BNOTP=0. 1150 

A1=AN0T/ CBCl-ACl > 

AA1=AN0TP/<BC1-AC1> 

B1=BN0T-A1*BC1 
BB 1 =BN0TP-AA 1 *BC 1 
A2=AN0T/ < BC2-AC2 ) 

AA2=AN0TP/<BC2-AC2) 

B2=BN0T-A2*BC2 

BB2=BN0TP-AA2*BC2 

1 WRITE<1, 1000) 

READ<1,*) ANS 
DO FOR 1=1,6 

IF<ANS.EQ. COMMANDED) THEN 

GOTO <100,200,300,400,500,600) I 

END IF 

ENDDO 

GOTO 1 

C BACKSPACE 

100 CONTINUE 

MTUX=MTUR ; ISTAT=MACTP< 15B , D , MTUX , D , D> 

GOTO 1 

C-.... REWIND 

200 CONTINUE 

MTUX=MTUR ; I STAT=MAGTP < 1 3B , D , MTUX , D , D ) 
IF<ISTAT.NE.O>THEN 
GOTO 999 
ENDIF 
GOTO 1 

READ---. 

00 DO FOR 1=1,1024 
IBUF1<I>=0 
ENDDO 

I DUM=0 i MTUX=MTUR j K0UNT=0 

2 CONTINUE 

PRINT’ <”^I>o You Want to READ the HEADER:”)’ 
READ’(A.l)* ,RAND 

IF<RAND.EQ. ’Y’ -OR.RAND.EQ. ’y’)then 
ISTAT=MAGTP< OB , I BUF2 , MTUR , 80 , I DUM) 



WRITE<1 ,1001>IDUI1 

WRITE< 1 , ’ <20A4) ’ ) < IBUF2< I ) , 1=1 , 40> 

END IF 
IJ=0 

4 CONTINUE 

IFCIBLK.EQ. 1000>G0 TO 17 
I STAT=MAGTP < OB , I BUF 1 , MTUR ,60,1 DUM > 

I F < ISTAT . EQ . 3>THEN 
PRINT’ <” EOF DETECTED”)’ 

GO TO 17 
ENDIF 

IBLK=IBLK+1 
DO 16 1=1,56,5 
IRE=IBUF1<I) 

IRE1=IBUF1<I+1> 

IRE2=IBUFl<I+2) 

IRE3=IBUFl(I+3> 

IRE4=IBUFl<I+4) 

RE1=FL0AT<IRE> 

RE2=FL0AT<IRE1) 

IF< IRE4.GT. 500- AND- IRE4.lt. 649)THEN 
JKI=1 

ELSEIF<:iRE4.GT.906.AND.IRE4.LT. 1052>THEN 
JKI=2 

ELSE 

GO TO 16 
ENDIF 

E1=CRE1>*TA1+TB1 

KEl=El+0.5 

E2=<RE2)*TA2+TB2 

KE2=E2+0.5 

EST= C F LOAT CIRE2)>*0. 032667-4 . 1 87703 
DEL=<FL0AT<IRE3)>*0. 02276 1 -0 . 425054 
ETOT=DEL+EST 
IRE5=ET0T*10 

IF<IRE5.LT.l .0R.IRE5.GT.750>IRE5=1 

ITAc(IRE5>=ITAC(IRE5)+1 

EX=59.7-ET0T-5.8 

IF<ET0T.LT.30.5.0R.ET0T.GT.54.5)G0 TO 16 
IFCKE1.lt. 30. OR. KE1.GT.125>G0T0 16 
IFCKE2.LT.30.0R.KE2.GT.125)GOTO 16 
I J K= < E TOT-30 . 5 > /2 . 0+ 1 
IFCIJK.LT.1.0R.IJK.GT.14>G0 TO 16 
IF<IRE2.LT.150>G0 TO 161 
MI=CET0T**1 .72-EST** 1.72) +0.5 
IFCMI.lt. 296. OR. MI. GT.340)G0 TO 16 
IPARCMI)=IPARCMI>+1 
161 TEMP1=E1 
TEMP2=E2 

AM 1=238 . ♦E2/ C E 1 +E2 ) ; 

AM2=238.-AM1 
NIND=0 
GO TO 6 

E1=CA1+AA1*AM1)*CRE1)+B1+BB1*AM1 
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E2= < A2+AA2*AM2 ) * < RE2) +B2+BB2*AM2 
AM1=238 . ♦E2/ < E 1 +E2> 

AM2=238.-AM1 

6 Z 1 =92. *AM 1/238. 

Z2=92.-Z1 

IND=1 

TE1=EI 

CALL L0SS<E1,AM1,Z1,IND,E1F> 

TE2=E2 

IND=2 

CALL L0SSCE2,AM2,Z2,IND,E2F> 

ELL1=E1F-TE1 

ELL2=E2F-TE2 

E1=E1F 

E2=E2F 

AM1=238 - *E2F / C E 1 F+E2F ) 

AM2=238.-AM1 

7 CALL KINEM<;ET0T,E1F,E2F,AM1,AM2,ECM1,ECM2) 
ECM1=E1F 

ECM2=E2F 

CALL NEUT 

AM1=AM1-RNU1 

AM2=AM2-RNU2 

NIND=NIND+1 

IF<NIND.LT.2)G0 TO 5 

IFCNIND.LT.3>G0 TO 7 

9 CONTINUE 
AM1=RNU1+AM1 
AM2=RNU2+AM2 
Ml=AMl+0.5 
M2=AM2+0 . 5 
ENET,=E1N+E2N 
NET=ENET+ 0.5 
NE1=E1N + 0.5 
NE2=E2N +0.5 

IE1<JKI,IJK,NE1) = IE1<:JKI,IJK,NE1>+1 
IE2C JKI , IJK,NE2>=IE2<JKI , I JK,NE2)+1 
IF<NET.LT. 1 . OR. NET. GT. 250) GOTO 10 
NTKE<JKI , IJK,NET>=NTKE<JKI , IJK,NET)+1 

10 IF<M1.LT-1.0R.M1.GT.200)GOTO 11 
NMKJKI , I JK,M1>=NM1 <JKI , IJK,M1>+1 

11 IF<M2.LT.1.0R.H2.GT.200>G0T0 16 

NMl(aKI,IJK,M2>=NMl(JKI,IJK,M2)+l 

12 INE=<NET-120>/10 
IFCINE.LT.l.OR.INE.GT.8)G0TO 13 

13 IMl=<Ml-60)/2 +1 
IF(IM1.LT.1.0R.IM1.GT.50>G0T0 14 

IME< JKI , UK, IMl , INE>=IME<JKI , UK, IMl , INE)+1 
IMEA< JKI , IMl , INE>=IMEA< JKI , IMl , INE>+1 

14 IM2=CM2-60)/2 +1 
IF(IM2.LT.1.0R.IM2.GT.50>‘60T0 15 

IMECJKI , UK, IM2, INE>=IME(JKI , UK, IM2, INE) + 1 
IMEAC JKI , IM2, INE)=IMEA< JKI , IM2, INE)+1 

15 JM1^M2 



IFCM1.GT.M2>JM1=M1 
KMl=<JMl-115>/5+l 
IF<KM1.LT.1.0R.KM1.GT.11)G0T0 16 
IN2=<NET-125)/2 +1 
IF<IN2.LT.1.0R.IN2.GT.40>G0T0 16 
IMRE<JKI,IJK,IN2,KM1) = IMRE<:JKI,IJK,IN2,KM1> + 1 

16 CONTINUE 
J=J+1 
GOTO 4 

500 PRINT’C”* Want to SKIP File:*’)’ 

REAP ’<A1>’, RAND 

IF<RAND.EQ. ’Y’ -OR-RAND-EQ- ’y’)THEN 
PRINT’ <”*How Many Files to Skip:”)’ 

READ*,IFILENO 
DO 502 IAS=1 ,IFILENO 
ISTAT=MAGTP<:iOB, IBUFl ,MTUR,60, IDUM) 

502 CONTINUE 
END IF 

IFCRAND.EQ. ’N’ -OR-RAND-EQ. ’n ’ )THEN 
PRINT’ <’ ’*How Many BLOCKS to Skip:”)’ 

READ*,IBLKN 
DO FOR IAS=1,IBLKN 

ISTAT=MAGTP<16B, IBUFl ,MTUR,60 , IDUM) 

ENDDO 
ENDIF 
GO TO 1 

17 WRITE<1,*) 

WRITE <1, 1002) IBLK 
WRITEd,*) 

PRINT’ <”4iDo You Want to CONTINUE -Y:Yes;N:No: ” ) ’ 

READ*<A1)’,RAND 

IBLK=0 

IF<RAND.EQ. ’Y’ -OR- RAND-EQ- ’y’ )THEN 
GO TO 2 
ENDIF 
GO TO 1 

600 CONTINUE 

C 

C WRITE<9,1005)<IPAR<I),I=1»2048) 

DO FOR 

c 

C WRITE<6, 1004)K 

WRITEC6,*)’ SINGLE FRAGMENT KINETIC ENERGY OF DCTECTOR-1 
WRITE(6,1005)<IE1<1,K,I>.I=1>200) 

DO FOR 1*1,200 

lEKl ,K,I> = IE1C1,K, I)-IE1<2,K,I) 

IFCIEKl ,K,I).LT.0)IE1<1,K,I)=0 
ENDDO 

WRITEC7,*)’ SINGLE FRAGMENT KINETIC ENERGY OF DETECTOR-1 
WRITE<7,1005)<IE1<1,K,I),I=1,200) 

C WRITE(6,1003) 

WRITE<6,*>’ SINGLE FRAGMENT KINETIC ENERGY OF DETECTOR-2 
WRITE C6, 1005) <IE2C1,K, I), 1=1, 200) 

DO FOR 1=1,200 



IE2<1 ,K, I>=IE2(1,K, I>-IE2<2,K, I) 

IF<IE2<1,K,I)-LT.0)IE2<1,K,I>=0 

ENDDO 

WRITE<7,*)’ SINGLE FRAGMENT KINETIC ENERGY OF DETECTOR-2’ 
WRITEX7, 1005> < IE2< 1 ,K, I ) , 1=1 ,200> 

C WRITE C6, 1003) 

WRITE<6,*>’ TOTAL KINETIC ENERGY’ 

WRITE<6,1005)<NTKE(1,K, I >,1=1,250) 

DO FOR 1=1,250 

NTKE<1 ,K,I)=NTKE<1,K,I)-NTKE<2,K,I) 

IF<NTKE<1,K,I).LT.0)NTKE<1,K,1)=0 

ENDDO 

WRITE<7,*)* TOTAL KINETIC ENERGY’ 

WRITE<7, 1005><NTKE<1 ,K, I) , 1=1 ,250) 

SUM1=0.0 

SUM2=0.0 

SUM3=0.0 

DO FOR 1=150,220 
SUM1=SUM1+FL0AT<NTKE<1,K,I)*I) 

SUM2=SUM2+FL0AT<NTKE < 1 , K , I ) ) 

ENDDO 

RMEAN=SUM1/SUM2 
DO FOR 1=150,220 

SUM3=SUM3+NTKE <1,K,I)*<< RME AN-F LOAT < I ) ) ** 2 ) 

ENDDO 

SIGM=SUM3/SUM2 

SIGM=SQRT<SIGM) 

WRITE<10,1004)K 

WRITE<10,1009)NTKE<1 ,K,NINT<RMEAN)),SUM2,RMEAN,,SIGM 
C WRITE<6,1003) 

WRITECG,*)’ MASS DISTRIBUTION’ 

WRITE<6,1005)(NM1<1,K,I),I=1,200) 

DO FOR 1=1,200 

NMKl ,K,I)=NM1<1,K,I)-NM1<2,K,I) 

IF(NM1<1,K,I).LT-0)NM1<1,K,I)=0 

ENDDO 

WRITE(7,*)’ MASS DISTRIBUTION’ 

WRITE<7 , 1005) <NM1 ( 1 ,K, I ) , 1=1 ,200) 

ENDDO 

WRITE<8,*)’ MASS DISTRIBUTION AS A FUNCTION OF TKE AVERAG- EX ENE 
DO FOR K=l,8 

WRITEC8,*) ’KINETIC ENERGY VALUE=’ , 120+10*K 
DO FOR J=l,50 

IMEA<1,J,K)=IMEA<1,J,K)-IMEA<2,J,K) 

IF<IMEA<1 , J,K)-LT.0)IMEA(1,J,K)=0 
ENDDO 

WRITECS, 1005)<IMEA<1,J,K),J=1,50) 

ENDDO 

c GO TO 1111 

DO FOR IJ=1,14 
WRITE<6,1004)IJ 
DO 20 1=1 ,8 
SM=0.0 
SAM=0.0 



SIGH=0.0 
NER=I*10+i20 
WRITE<6, 1006>NER 
DO FOR K=l,50 

IMEC1,IJ,K,I>=IME<1,IJ,K,I)- IME<2,IJ,K,I> 

IF<IME(1,IJ,K,I>.LT.6>IME<1,IJ,K,I>=0 

ENDDO 

WRITE<6,1003> 

WRITE<6,1005><IME<1,IJ,K,I>,K=1,50> 

RM=60.5 

DO 18 J=l,50 

SM=SM+IME<1,IJ,J,I> 

SAM=SAM+RM* IME<1,IJ,J,I) 

RM=RM+2.0 

18 CONTINUE 
IFCSM.EQ.O.O>GOTO 20 
SAM=SAM/SM 
RJH=60.5 

DO 19 JJ=1,50 

S I GM=S I GM+ IME<1,IJ,JJ,I)*< R JM-SAM > **2 
RJM=RJM+2.0 

19 CONTINUE 

S I GM=SQRT < S I GM /SM ) 

WRITE<6, 1007)NER,SM, SAM, SIGN 

20 CONTINUE 
ENDDO 

DO FOR IJ=1,14 

WRITEC7,1004>IJ 

DO 23 1=1 ,8 

SM=0.0 

SAM=0.0 

SIGM=0.0 

NER=I*10+120 

WRITE<7, 1006)NER 

WRITEC7,1003> 

WRITE <7, 1005)<IME<2,IJ,K,I>,K=1,50> 

RM=60.5 
DO 21 J=l,50 
SM=SM+IME<2, IJ,J,I> 

S AM=S AM+RM* IME<2,IJ,J,I) 

RM=RM+2.0 

21 CONTINUE 
IF<SM.EQ.O.O>GOTO 23 
SAM=SAM/SM 
RJM=fe0.5 

DO 22 JJ=1,50 

S I GM=S I GM+ IME<2,IJ,JJ,I)*( R JM-SAM > **2 
RJM=RJM+2.0 

22 CONTINUE 

S I GM=SQRT < S I GM/ SM ) 

WRITE(7, 1007)NER,SM,SAM,SIGM 

23 CONTINUE 
ENDDO 

DO FOR LM=1,14 



WRITE<6,1004>LM 

DO 26 1=1,14 

SE=0.0 

SFM=0.0 

SIGE=0.0 

RI=I 

MASS=120+I*5 
WRITE<6,1008)MASS 
WRITE <6, 1003) 

DO FOR K=l,40 

IMRE<1,LM,K,I)=IMRE<1,LM,K,I)-IMRE<2,LM,K,I> 
IF(IMRE<l,LM,K,I).LT.O>IMRE<:i ,LM,K, I>=0 
ENDDO 

WRITE<6,1005><IhRECl,LM,K,I>,K=l,40> 

RE=120.5 
DO 24 J=l,40 
SE=SE+IMRE<1 ,LM, J,I> 
SFM=SFM+RE*IMRE(1,LM,J,I) 

RE=RE+2.0 

24 CONTINUE 
IFCSE.EQ.O.O>GOTO 26 
SFM=SFM/SE 
RJE=120.5 

DO 25 JJ=1,40 

SIGE=SIGE+IMRE< 1 , LM, J J , I >*<RJE-SFM)**2 
RJE=120.5 

25 CONTINUE 
SIGE=SQRT<SIGE/SE) 
WRITE<6,1009>MASS,SE,SFM,SIGE 

26 CONTINUE 
ENDDO 

DO FOR LM=1,14 

WRITE<7,1004>LM 

DO 29 1=1,14 

SE=0.0 

SFri=0.0 

SIGE=0.0 

RI=I 

MASS=120+I*5 

WRITE<7,1008>MASS 

WRITE<7,1003) 

WRITEC7, 1005)<IMREC2,LM,K,I) ,K=1 ,40) 

RE=120.5 

DO 27 J=l,40 

SE=SE+IMRE<2,LM,J,I) 

SFH=SFM+RE*IMRE<2,LM, J, I> 

RE=RE+2-0 

27 CONTINUE 
IF<SE.EQ.O.O>GOTO 29 
SFM=SFM/SE 
RJE*120.5 

DO 28 JJ=1,40 

SIGE=SIGE+IMRE<2,LM,JJ,I)*<RJE-SFM>»*2 

RJE=RJE+2.0 



28 CONTINUE 
SIGE=SQRT(SIGE/SE) 

WRI TE < 7 , 1 009> MASS , SE , SFM , S I GE 

29 CONTINUE 
ENDDO 
GOTO 1 


C-..ERR 

999 WRITE<1,*> »♦****» 

WRITE Cl,*) ’ERRNO=» ,ISTAT 
GOTO 1111 

C... EXIT 

400 GOTO 1111 
2345 CONTINUE 
1011 ICNUM=ICNUM+1 

IFCISTAT.EQ.3>THEN 
WRITECl,*) ’ END OF FILE’ 

WRITECl,*) ’ BLOCK=’ ,IBLK, ’ I STAT= ’ , I STAT 
ICNUM=ICNUM+1 

IF<ICNUM.LT.INUM)GOTO 2345 

GOTO 1 

ELSE 

GOTO 999 
ENDIF 

1111 CLOSECG) 

CL0SEC7> 

CL0SEC3> 

CL0SE<9) 

CLOSECIO) 

STOP 

END 

SUBROUT INE KINEMC ETOT , E 1 , E2 , AM 1 , AM2 , ECM 1 , E CM2 > 

COMMON/AA/EP 

COMMON/EE/THETAl ,THETA2 

AM=238.0 

AP=4.0 

PI=4.0*ATANC1.0) 

VP=i .389*SQRTCEP/AP> 

ETAR=59-7-ET0T 
C VCN=AP*VP/AM 

VCN=1 .389*SQRT<ETAR/AM> 

VLl=l .389*SQRT<E1/AM1) 

VL2==1 .389*SQRT<E2/AM2) 

ANG1=AC0SCVCN/VL1 ) 

ANG2=AC0S C VCN/ VL 2 ) 

VCM 1 = VL 1 **2-VCN**2 

VCM2=V/L2**2-VCN**2 

ECM 1 =AM 1 *VCM 1 / 1 . 929321 

ECM2=AM2*VCM2/1 .929321 

RETURN 

END 

SUBROUTI NE LOSS < ERES , AMF , ZF , IND , E F > 



C*********THIS SUBROUTINE CALCULATES THE ENERGY LOSS SUFFERED***** 
C*********«***«BY FISSION FRAGMENTS***************** 

COMMON/ BB/SU < 20 > 

COMMON/CC/SC < 20 > 

COMMON/EE/THETAI ,THETA2 

COMMON/RR/BACK , THICK 

EL=0-0 

EF=0.0 

NCN=0 

PI=d.0*ATAN<1.0> 

NCN2=0 

GOTO <20,24>IND 
20 THETA=THETAl-45. 

THETA=THETA*P I / 1 80 . 

TU= . 05/C0S(THETA) 

GOTO 69 

24 THETA=THETA2-45. 

THETA=THETA*PI/180. 

TU= . 05/C0S< THETA) 

69 EBM=ERES/AMF 

I=EBM/0.1 
ZALF=2. 

S1=SU<I>-<SU<I)-SU<I+1))*<EBM-0.1*I)/0.1 
63 V=1.389*SQRT<ERES/AMF) 

V0=0.2188 

VR1=V/(V0*ZALF**.65) 

VR2=V/ ( VO*ZF** - 65) 

AZ1=1 .035-0. 4*EXP<-0. 16*2. ) 

AZ2=1 .035-0. 4*EXPC-0.16*ZF) 

GAMA 1=1. 0-AZ 1 *E X P C -0 . 879*UR1 ) 

GAMA2=1 .0-AZ2*EXP(-0.879*VR2) 

GAMA 1=< GAMA 1*ZALF)**2 | 

GAMA2= < GAMA2*Z F ) **2 

S=S 1 *G AM A2 / GAMA 1 

EL=S*TU 

EF=ERES+EL 

IF(IND.EQ.1)G0T0 74 j 

IF(IND.EQ.3)GOT0 74 
72 ERES=EF 

TU= . 0 1 /ABS < COS < THETA) ) 

IND=IND+1 

75 EBM=ERES/AMF 

I=EBM/0.1 

S1=SC<I)-<SCCI)-SC<I+1))*<EBM-.1*I)/0.1 
GOTO 63 
74 RETURN 

END 

SUBROUTINE NEUT 
COMMON/DD/EX 

*»*****************************************************************1 
C THIS PROGRAM CALCULATES THE NUMBER OF NEUTRONS EMITTED FROM 

C NASCENT FISSION FRAGMENTS AT LOW EXCITATION ENERGY OF TARGET 

C THE DATA USED IS FROM THE PAPER OF BISHOP ETAL . 

*#*********4HHH(.************************-»******-»****-»*»*************^ 



COMMON/NEUT/ECMl ,ECM2,E1N,E2N,RNU1 ,RNIJ2,AM1 ,AM2 
DIMENSION RNU(2) 

AMI = ECM2*238.0/<:ECM1+ECM2> 

AM2 = 238.0-AMl 
IFCAM1.GT.AM2>THEN 
MASS = AMI + 0-5 
LIND=1 
ELSE . 

MASS=AM2 
LIND=2 
END IF 

SLOPE 1=0.0517885 
SL0PE2=0. 0542135 
1=1 

IF<MASS.LE.120>THEN 

IF(EX.GT.5.7>THEN 

RINT = -3.3459896 + O. 0497349* < EX-5 -7) 
ELSEIFCEX.LT.5.7>fHEN 

RINT = -3-3459896 - O- 0497349* <EX-5. 7 > 

ELSE 

RINT = -3.3459896 
ENDIF 

RNU<I) = SL0PE1*MASS + RINT 

ENE = <FLOAT< MASS-50 )> * 0-004117647 + 1-1268382 
IF(LIND.EQ.1>E1N=ECM1+ENE*RNU<I> 

I F < L I ND . EQ . 2 > E2N=ECM2+ENE*RNU C I > 

1=1+1 

ENDIF 

I F ( MASS . GT . 1 32 ) THEN 
IF<EX.GT.5.7)THEN 

RINT = -5.350377 + 0-0514494*<EX-5. 7> 
ELSEIF<EX-LT.5.7>THEN 

RINT = -5.350377 - 0-0514494*<EX-5.7> 

ELSE 

RINT = -5.350377 
ENDIF 

RNUCI) = SL0PE2*MASS + RII^T 

ENE = <FLOAT<MASS-50>> * 0-004117647 + 1-1268382 

IF<LIND.EQ.1)E1N=ECM1+ENE*RNUCI) 

IF<LIND.EQ.2>E2N=ECM2+ENE*RNU(:i> 

1=1 + 1 
ENDIF 

I F < M ASS - GT . 1 20 . AND • MASS . L E - 1 32 > THEN 
IF<EX.GT.5.7>THEN 

RINTl = -3-3459896 + 0. 0475725* (EX-5.7> 
ELSEIF<EX.LT.5.7)THEN 

RINTl = -3.3459896 - 0 - 0475725*<EX-5 - 7> 

ELSE 

RINTl = -3.3459896 
ENDIF 

IF<EX.GT.5.7>THEN 

RINT2 = -5.350377 + 0.0492125*<EX-5.7> 
ELSEIF(EX.LT.5.7>THEN 

RINT2 = -5.350377 -0.0492125*<EX-5.7> 
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ELSE 


RINT2 = -5.350377 
ENDIF 

RPl = SLOPEl * 120-0 + RINTl 
RP2 = SL0PE2 * 132.0 + RINT2 
RINT3 = C RPl *11-0 - RP2 * 10-0) 

RNU<I> = CCRP2-RP1>/12.0>*MASS + RINT3 

ENE = <FL0AT<MASS-50>> * 0-004117647 + 1-1263382 

IFCLIND-EQ. 1>E1N=ECM1+ENE*RNUCI> 

I F < L IND - EQ - 2 > E2N=:ECM2+ENE*RNU< I ) 

1 = 1+1 

ENDIF 

IF<I.GT.2)G0 TO 2 
IFCLIND-EQ-DTHEN 

MASS = AM2 + 0-5 
LIND=2 

ELSE 

MASS = AMI + 0-5 
LIND=1 
ENDIF 
GO TO 1 

2 IF<LIND.EQ.1>THEN 

RNU1=RNU<2) 

RNU2=RNUC1> 

AM2=238 - 0*E 1 N/ < E 1 N+E2N > 

AM1=233.0-AM2 

ELSE 

RNU1=RNU<1> 

RNU2=RNU<2> 

AM1=238.0*E2N/<E1N+E2N) 

AM2=238.0-AM1 
ENDIF 

EPM1=E1N/AM1 
EPM2=E2N/AM2 
E1N=E1N+RNU1*EPM1 
E2N=E2N+RNU2*EPM2 
RETURN 
END 
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APPENDIX 3 


THIS PROGRAM CALCULATES THE NUMBER OF NEUTRONS EMITTED FROM 
COMPOUND NUCLEUS TO THE SADDLE POINT 
PROGRAM PRE SADDLE 

DIMENSION EXCN<20> ,GAMMAF<20) ,GAMANC20> ,GAMAGC20> ,GAMAT<20) , 

- ICHANC<20),BS<45),BC<45) 

COMMON/ PH I / P I , P I SQ 
COMMON/MCS/A,Z,S 

C0MM0N/C0NS/DEF<45) ,ELDM<45> ,DELS<45> ,EDEF<45) , ALR0<45) ,RH0<45) 
COMMON/CONSl /AS, CAPA , Cl , W , ALDM , ESO , ECO 
COMMON/BS/ < 45 > , BC < 45 ) 

COMMON/MOM/XJPAR,XJPER 
COMMON/ A/ A 1 , A2 , A3 

DATA BC/O . 9976 IDO , O . 99849D0 , 0 . 999 1 6D0 , 0 - 99963D0 , O - 99991 DO , 

1 . OOOOODO , O . 99991 DO , O . 99965DO , O - 9992 1 DO , O - 99862DO , 

O . 99786DO , O - 99695DO , O - 99589DO , O . 99467D0 , 0 . 99330DO , 

0.991 78D0 , 0 . 990 1 1 DO , O . 98330D0 , O . 98634D0 , 0 . 98423D0 , 

O . 98 1 97D0 , O . 97956D0 , O . 9770 1 DO , O . 97430D0 , 0 . 97 1 44DO , 

O . 96843D0 , O - 96526D0 ,0-961 93D0 , 0 . 95845D0 , 0 - 95480D0 , 

O . 95099D0 , O . 94702D0 , 0 . 94289D0 , 0 - 93859D0 , 0 . 934 1 2D0 , 

O . 92949D0 , O . 92469D0 , O . 9 1 973DO , O . 9 1 460D0 , O . 90932D0 , 

O . 90387D0 , 0 . 89826D0 , 0 . 89254D0 , 0 . 88666D0 , 0 . 88065D0/ 

DATA BS/1.00497DO,1.00312DO,1-00172DO,1.00075DO,1.00018DO, 

1 - OOOOODO , 1 . 000 1 8D0 , 1 . 00070DO , 1 . 00 1 54D0 , 1 . 00269DO , 

1 - 004 1 5D0 , 1 - 0058800 , 1 . 00789D0 ,1-0101 7D0 ,1-01 270D0 , 

1 - 0 1 549D0 , 1 - 0 1 85 1 DO , 1 . 02 1 78D0 , 1 - 02528D0 , 1 - 0290 1 DO , 

1 - 03296D0 , 1 - 03714D0 , 1 - 041 54D0 , 1 - 04616D0 , 1 - 05099D0 , 

1 - 05604D0 ,1-061 30D0 , 1 - 06678DO , 1 - 07247D0 , 1 - 07837D0 , 

1 - 08448D0 , 1 - 09080DO , 1 - 09732D0 ,1-1 0406D0 ,1-11 099DO , 

1 - 1 1813D0, 1 - 12546D0, 1 - 13299D0, 1 - 14071D0, 1 - 14861DO, 

1 - 15670D0, 1 - 16495D0, 1 - 17336D0, 1 - 18193D0, 1 - 19064D0, 

DATA X J PAR/ 1-11 578 , 1 - 09086 , 1 - 06688 , 1 - 04378 ,1-02151,1- 00000 , 


0 - 9792 1,0- 95908 , 0 - 93958 , 0 - 92067 , 0 - 90230 , 0 - 88445 , 
0.86708,0.85018,0.83370,0.81763,0.80196,0.78665, 
0.77170,0.75709,0.74281,0.72884,0.71519,0.70184, 
O . 68879 , 0 . 67604 , 0 . 66359 , 0 - 65144 , 0 - 65960 , 0 . 62806 , 
0.61685,0.60596,0.59541,0.58523,0.57541,0.56599, 
O . 55698 , 0 . 54842 , O . 54032 , O - 53272 , 0 . 52564 ,0.51914, 
0.51323,0.50797,0.50341/ 

DATA X JPER/0 . 95952 , 0 . 96552 , 0 . 97272 , 0 . 98083 , O - 98992 , 1 - OOOOO , 

1.01107.1. 02312 . 1 . 03618 . 1 - 05025 , 1 .06535 , 1 . 08151 , 

1 . 09875 .1.1 1709 . 1 . 13658 . 1 - 15726 ,1-17916,1. 20233 , 

1 . 22683 . 1 . 25272 . 1 . 28005 . 1 - 30890, 1 . 33933 ,1.37 143 , 
1.40527,1.44095,1.47856,1.51820,1.55997,1.60400, 

1 . 65039 .1 . 69927 . 1 . 75078 . 1 - 80506 , 1 - 86225 , 1 - 92251 , 

1.98599.2.05288.2.12334.2- 19757,2.27576,2.35812, 



2.44485,2.53618,2.63235/ 

PRINT’ <”*PLEASE GIVE THE BINDING ENERGY OF THE NEUTRONS 
RE AD*, BN 

PI=<ATAN(1 .0>>*4.0 
PISQ=PI*PI 

PRINT’ C”*PLEASE GIVE THE EX- ENERGY OF THE NUCLEUS:”)’ 

READ*, EX 

ALD=3.293D0 

HRED=6.588D-22 

IN=1 

PRINT’ < ”«PLEASE GIVE THE RANDOM NUMBER SEED HERE:”)’ 

READ*, IX 

IHIS=1 

DT=1.0D-21 

T=0.6 

EXCN=EXCN-2.0 

GAMAF=GAMAFCI) 

AS=19.008 

CAPA=2.S4 

Cl=0.731 

W=0 . 185*A** < 1 . 0 /3 . 0 ) 

ALDM=A*0. 176*<1 .0-A**< 1 .0/3-0) ) 
WRITE<2,1000)A,Z,S,EX,ALDM 

ESO=AS* < 1 . 0-CAPA*< ( A-2 . 0*Z ) / A ) **2) *A** < 2 - 0/3 . 0) 
EC0=C1*Z*Z/A**<1. 0/3.0) 

DO 3 1=1,45 
DEF<I)=<I-1)*0.02+0.9 

ELDM(I)=ESO*<BS<I>-1.0)+ECO*<BC<I)-l.O) 

THETA=<DEFCI)-1 .0)*A**C1 .0/3.0)/0.24 

DELS<I)=S*SIN<THETA)/THETA 

EXDEF<I)=ELDM<I)+DELS(I) 

EXCI)=EXCN-EDEF<I) 

WRITE<2, 1010)1 ,DEF<I),ELDM<I),DELS<I),EDEF<I),EX<I) 
CONTINUE 

GAMN=CALL TAUN<EXE ,GAMAN) 

GAMF=CALL TAUF<EXE ,GAMAF) 

GAMG=CALL TAUG<EXE ,GAMAG) 

GAMT=GAMN+GAMF+GAMG 

ICH=1 

PSRAND=RNDM( I X , IN) 

TEMP=CEX/ALD>**0.5 

BRF=GAMF/GAMT 

BRN=GAMN/GAMT 

PROBF=BRF*GAMT*DT/HRED 

PROBN=BRN*GAMT*DT/HRED 

IF<PR0BF-PSRAND)150,150,100 

WR ITE* , T , EX ,6AMT , PROBF , PROBN , PSRAND 

ICHANC< ICH)=ICHANC< ICH)+1 

PRINT*, IHIS 

GO TO 500 

I F (PSRAND. GT. PROBF. AND. (PSRAND-LT.CPROBF+PROBN)) GO TO 200 

IFCPSRAND.GT. (PROBF+PROBN))GO TO 300 

WR I TE <*, 1 030 ) T , E X ,GAMT , PROBF , PROBN , PSRAND 

FORMATC lOH, FISSION, /,6D11 .4) 



1030 FORMATClOH, NEUTRON, /,6D11.4> 

EX=EX-BN-2.0*TEMP 
IF(EX.LE.10.0D0>G0 TO 700 
T=T+DT 
ICH=ICH+1 

GAMF=GAMAF C I+2*< ICH-1 ) > 

G AMN=G AM AN < I +2* < I CH- 1 > > 

G AMG=GAM AG < I +2* < I CH- 1 ) > 

G AMT=G AMF +G AMG+G AMN 
GO TO 60 

300 T=T+DT 

GO TO 60 

500 IHIS=IHIS+1 

IF (IHIS-lOO) 10,600,600 

600 PRINT 1040, (ICHANC< I), 1=1,10) 

IF<IHIS.EQ. 100>HIS=IHIS 

1040 F0RMAT<10I8> 

NEUTR=0 
DO 680 1=1,10 

680 NEUTR=NEUTR+ICHANC<I)*I 

ANEUTR=NEUTR 
RNEUTR=ANEUTR/HIS 
WRITE<*, 1050>RNEUTR 

1050 FORMAT <26H, AVERAGE CHANCE NEUTRON IS =,F10-2) 

STOP 

700 ICHANC<;i0) = ICHANC<10) + l 

GO TO 10 
END 

FUNCTION RNDM<IX,IN) 

INTEGER*4 JM,MM,IX,IC 
IF(IN>30,4O,30 

30 IN=0 

JM=2**10+3 

MM=2**(20> 

FM=MM 

40 IC=IX»JM/MM 

IX=IX*JM-IC*MM 

FX=IX 

XD=FX/FM 

RNDM=ABS(XD) 

RETURN 

END 

SUBROUTINE TAUF<EXE ,6AMAF ,GAMAT) 

IMPLICIT DOUBLE PRESCISION <A-H,0-Z) 

COMMON/PHI/PI ,PISQ 
COMMON/MCS/A Z S 

C0MM0N/C0NS/DEF<:45) ,ELDM<45) ,DELSC45) ,EDEF<45) , ALRO<:45> ,RH0<:45) 
COMMON/CONSl / AS , CAPA , Cl , W , ALDM , ESO , ECO 
C0MM0N/BS/(45) ,BC(45) 

COMMON/A/Al ,A2,A3 
DIMENSION EX(45> 

DO 10 1=1,45 

SPIN=0 

DLS=DELS(I> 



EX1=EX(I> 

CALL LEVDEN<A,ALDM,DLS,W,EX1 ,SPIN,WE,SIGSQ,RH0E,RH01 ,ALR01> 
ALR0<I>=ALR01 

WRITE<2, 1000)1 ,RHO<I> ,ALRO<I> 

F0RMATCI3,E13.5,E13.5) 

CONTINUE 

IIMAX=43 

IGS=6 

ISP=6 

DO 20 11=2, UMAX 

IF<:RHOCII+l).LT.RHO<;iI).AND.RHOai-l>-LT.RHO<II)>GO TO 21 
GO TO 22 

IF<RHO< I I ) .GT .RHO< IGS) ) IGS=I I 
CONTINUE 

IF<RH0<II+1).GT.RH0<II).AND.RH0CII-1> -GT.RHO<;iI)>GO TO 23 
GO TO 20 

IFCISP-6)24,24,25 

ISP=II 

GO TO 20 

IF<RHO<II) .LT.RHO<ISP)>ISP=II 
CONTINUE 

WRITE<2,*>’IGS=’ ,IGS 
WRITE<2,*> ’ISP=’ ,ISP 
DLS=DELS<ISP) 

EX1=EX<ISP) 

EX2=EX1-1.0 
EX3=EX2-1.0 
EX4=EX3-1.0 
EX5=EX4-1 .0 
EX6=EX5-l-0 

CALL LEVDEN<A,ALDM,DLS,W,EX1,SPIN,WE,SIGSQ,RH0E,RH01,ALR01) 

CALL LEVDEN(A,ALDM,DLS,W,EX2,SPIN,WE,SI GSQ , RHOE , RH02 , ALR02> 

CALL LEV/DEN<A,ALDM,DLS,W,EX3,SPIN,WE,SIGSQ,RHOE,RH03,ALR03) 

CALL LEV/DEN<A, ALDM ,DLS,W,EX4,SPIN,WE,SI GSQ , RHOE , RH04 , ALR04 > 

CALL LEVDEN<A, ALDM, DLS,W,EX5, SPIN, WE, SIGSQ, RHOE, RH05,ALR05) 

CALL LEVDEN<A, ALDM, DLS,W,EX6, SPIN, WE, SIGSQ, RHOE, RHO6.,ALR0fe> 
ALROO=ALRO<IGS) 

GAMA=EXP<ALR01-ALR00)/2.0 

GAMA=GAMA+EXP<ALR02-ALR00) 

GAMA=GAMA+EXP<ALR03-ALR0O) 

GAMA=GAMA+EXP<ALR04-ALR00> 

GAMA=GAMA+E X P ( ALROS-ALROO) 

GAM A=GAMA+E X P < AL R06- ALROO ) 

GAMA=GAMA/2./PI 
GAMAF=GAMA 
GAMT=GAMAF*3.2E21 
TAUF='l .O/GAMT 
WRITEC2, 1010>TAUF 
WRITE<2, 1020)GAMAF,GAMT 

FORMATCTH GAMAF=,E13.5,3HMEV,5X ,E13.5,3H/S> 

F0RMAT<7H TAUF=,E13.5,3H/S) 

RETURN 

END 

SUBROUT I NE L E VDEN < A , ALDM , DELS , W , U , SP I N , WE , S I GSQ , RHOE , RHOE I , ALRH > 



COMMON/PHI /PI ,PISQ 
COMMON/ A/A 1 ,A2, A3 

IMPLICIT DOUBLE PRESCISON<A-H,0-Z) 

WE=1.0 

SIGSQ=0.0 

RH0E=1.0 

RH0EI=1.0 

ALRH=0.0 

IF<lJ.LE.O.O>RETURN 

TP=SQRT<U/ALDM) 

10 T=TP 

T1=T 

T2=T+0.005 

U1=ALDM*T1*T1-DELS 

U2=ALDM*T2*T2-DELS 

XX1=PI»W*T1 

IF<XX1 .GT.15.0)G0 TO 700 
XX2=PI*W*T2 

U1=U1+DELS*C<XX1>**2*C0SH<XX1)/CSINH<XX1))**2) 

U2=U2+DELS*C<XX2>**2*C0SH<XX2>/<SINH<XX2))**2) 

700 CONTINUE 

TP=T1+(U-U1 )*<T2-T1 >/<U2-Ul > 

IFCABS<T-TP) .GT-0.005)GO TO 10 
T=TP , 

S=2.0*ALDM*T 

XXX=PI*W*T 

IF<XXX.GT.15-0>G0 TO 99 

S=S+DELS/T*<<XXX>**2*C0SH<XXX)/<SINH<XXX>)**2-XXX/SINH(XXX)> i 
99 AEFF=S/2.0/T 

GEFF=6.0*AEFF/PISQ 

EM2BAR=0.240*A**<2.0/3.0) I 

SISSQ=EM2BAR*GEFF*T 1 

C=SQRT<PI>/12-0/AEFF**0.25/U**C5.0/4.0> 

Cl=l .0/(2.0*PI*SIGSQ)**0.5 
ALRH=S+DL0GCC>+DL0G<C1 ) 

WE=EXP(S)»C 

RH0E=WE*C1 

RH0EI=RH0E*C2.0*SPIN+l>*EXP<-SPIN*<SPIN+l>/2.0/SI6SQ) 

RH0EI=RH0EI/2.0/SIGSQ 

RETURN 

END 

SUBROUTINE TAUNCEXE.GAMAN) 

IMPLICIT DOUBLE PRECISI0N<A-H,0-Z) 

COMKON/PHI/PI ,PISQ 
COMMON/MCS/A , Z ,S 

COMMON/CONS/DEFC45),ELDM<45>,DELS<45),EDEF<45) ,ALR0<:45>,RH0(45> 
COMMON/CONSl / AS , CAPA , Cl , W , ALDM , ESO , ECO 
COMMON/BS/ C 45 ) , BC C 45 ) 

COMMON/MOM/XJPAR,XJPER 

COMMON/A/Al ,A2, A3 

Bl=31.05 

B2=-25.91 

B3=342.4 

B4=21.89 



B5=0.223 
B6.=0 - 673 
B7=617.4 
AMAS=A 

A1=B1*AMAS**<-1 .0/3.0>+B2 

A2=B3* AMAS** C-l-0/3.0> +B4»AMAS** C 2 - 0/3 - 0 ) 

A3=B5* AMAS** < - 1 . 0/3 . 0 ) +B6*AMAS** < 2 . 0/3 - 0 ) +B7 

DO 10 1=1,45 

SP.IN=0 

DLS=DELSCI) 

EX1=EXCI> 

CALL LEVDEN<A,ALDM,DLS,W,EX1 ,SPIN,WE,SIGSQ,RHOE,RH01 ,ALR01> 
ALRO<I)=ALR01 

WRITE<:2,10O0>I,RHO<I),ALR0<I) 

1000 F0RMAT<I3,E13.5,E13.5> 

10 CONTINUE 

IIMAX=43 
IGS=6 
ISP=6 

DO 20 11=2, UMAX 

IF(:RHO<II+l).LT.RHO<II).AND.RHO<:iI-l).LT.RHO<ir>>GO TO 21 
GO TO 22 

21 IF<RHO<II) .GT.RHO<IGS))IGS=II 

22 CONTINUE 

IF<:RH0(II + 1).GT.RH0<II).AND.RH0<II-1>.GT.RH0<II>>G0 to 23 
GO TO 20 

23 IF<ISP-6)24,24,25 

24 ISP=II 
GO TO 20 

25 IF<RHO(II>.LT.RHO<:iSP>>ISP=II 

20 CONTINUE 

DO 30 11=1,45 
DLS=DELS<II> 

EX1=EX<II) 

EX2=EXl-BN-0.5 

EX3=EX2-1.0 

EX4=EX3-1.0 

EX5=EX4-1.0 

EX6=EX5-1.0 

CALL LEVDEN<A,ALDM,DLS,W,EX1,SPIN,WE,SIGSQ,RH0E,RH01,ALR01) 
CALL LEVDEN(A,ALDM,DLS,W,EX2,SPIN,WE,SIGSQ,RH0E,RH02,ALR02) 
CALL LEVDEN<A,ALDM,DLS,W,EX3,SPIN,WE,SIGSQ,RH0E,RH03,ALR03) 
CALL LE\/DEN<A,ALDM,DLS,W,EX4,SPIN,WE,SIGSQ,RH0E,RH04,ALR04) 
CALL LEVDEN<A,ALDM,DLS,W,EX5,SPIN,WE,SIGSQ,RH0E,RH05,ALR05> 
CALL LEVDEN<A,ALDM,DLS,W,EX6,SPIN,WE,SIGSQ,RH0E,RH06,ALR06> 
SIG1=SIGR<0.5D0) 

G A M A=D E X P < A L R02- A L RO 1 ) *0 - 5*93 1 - 0*S I G 1 / 1 0 - 0 
G AM A=GAM A+DE X P < ALR03-ALR0 1 > *1 . 5»931 . 0*SI GR < 1 - 5D0 ) / 1 0 - 0 
GAM A=G AM A+DE X P < AL R04- ALRO 1 ) *2 . 5*93 1 - 0*S I GR < 2 - 5D0 > / 1 0 - 0 
G AM A=G AM A+DE X P < ALR05-ALR0 1 ) *3 . 5*93 1 - 0*S I GR < 3 - 5D0 > / 1 0 - 0 
GAMA=GAMA+DEXP< ALR06-ALR01 )*4 - 5*931 - 0*SIGR<4 . 5D0) /lO. 0 
G AM < 1 1 ) =GAM A/ < P I + 1 97 . 0 ) **2 
30 CONTINUE 

GAMAN=0.0 



PR0B=0 - 0 

DO 31 11=1 , ISP 

G AH AN=G AMAN+GAM C I I ) *RHO < I I ) 

PROB=PROB+RHO < 1 1 > 

31 CONTINUE 

GAHAN=GAMAN/ PROS 
GAMT=GAMAN*3 . 2E21 
TAUN=1 -O/GAMT 
RETURN 
END 

FUNCTION SIGR 

CC****»***NEUTRON CROSS SECTIONS- CHATTERJEE ETAL *********** 
CC******-k-******PRAMANA 16< 1981 )391***********-»***** 

COMMON/A/Al ,A2, A3 

IMPLICIT DOUBLE PRECISI0N<:A-H,0-Z> 

S I GR=A 1 ♦EN+A2+A3/ EN 

RETURN 

END 

SUBROUTINE TAUG<EXE ,GAMA6) 

IMPLICIT DOUBLE PRECISION<A-H,0-Z> 

COMMON/PHI/PI ,PISQ 
COMMON/MCS/A,Z,S 

COMMON/CONS/DEFC45> ,ELDM<d5> ,DELS<45> , EDEF<45> , ALRO<45> ,RHO<453 
COMMON/CONS1/AS,CAPA,C1 ,W,ALDM,ESO,ECO 
C0MM0N/BS/<45> ,BC<45) 

COMMON/MOM/ XJ PAR, XJPER 
COMMON/A/Al ,A2,A3 

CC»********PARAMETERS FOR GAMMA INVERSE CROSS SECTION CALCULATIONS***** 
CC*********************PRAKASH ETAL- PRL 60(1988)1630************* 
GAM1=6.0 

ER=34.0/A**(1 .0/6.0) 

DO 10 1=1,45 

SPIN=0 

DLS=DELS(I) 

EX1=EX(I) 

CALL LEVDEN<A,ALDM,DLS,W,EX1,SPIN,WE,SIGSQ,RH0E,RH01,ALR01) 

ALR0(I)=ALR01 

RH0<;i)=RH01 

WRITE(2, 1000)1 ,RHO(I),ALRO(I) 

1000 F0RMAT<I3,E13.5,E13.5) 

lO CONTINUE 

IIMAX=43 
IGS=6 
ISP=6 

DO 20 11=2, UMAX 

IF(RH0<II+1).LT.RH0<II).AND-RH0(II-1).LT.RH0(II))G0 TO 21 
GO TO 22 I 

21 IF(RHO(II).GT.RHO(IGS))IGS=II I 

22 CONTINUE 

IF<RH0<II+1) .GT.RH0<II).AND.RH0(II-1).GT.RH0CII))G0 TO 23 
GO TO 20 

23 IF<:iSP-6)24,24,25 

24 ISP=II 
GO TO 20 



25 IF(RHO<II>.LT.RHO<ISP>>ISP=II 

20 CONTINUE 

DO 30 11=1,45 
DLS=DELS<II> 

EX1=EX<II> 

EX2=EXl-0.5 

EX3=EX2-1.0 

EX4=EX3-1.0 

EX5=EX4-1.0 

EX6=EX5-1.0 

CALL LEVDEN<A,ALDM,DLS,W,EX1,SPIN,WE,SIGSQ,RH0E,RH01,ALR01> 
CALL LEVDEN<A,ALDM,DLS,W,EX2,SPIN,WE,SIGSQ,RH0E,RH02,ALR02> 
CALL LEVDEN<A,ALDM,DLS,W,EX3,SPIN,WE,SIGSQ,RH0E,RH03,ALR03) 
CALL LE\/DEN<A,ALDM,DLS,W,EX4,SPIN,WE,SIGSQ,RH0E,RH04,ALR04> 
CALL LEVDEN<A,ALDM,DLS,W,EX5,SPIN,WE,SIGSQ,RH0E,RH05,ALR05> 
CALL LEVDEN<A,ALDM,DLS,W,EX6,SPIN,WE,SIGSQ,RH0E,RH06,ALR06) 
SIG1=SIGR<0.5D0> 

G AHA=D E X P ( A L R02- A L RO 1 ) *0 . 5**2»S I G I / 1 0 - 0 
G AM A=G AM A+DE X P ( AL R03- ALRO 1 > * 1 - 5**2*S I GR < 1 - 5D0 ) / 1 0 - 0 
GAMA=G AMA+DE X P < ALR04-ALR0 1 ) *2 . 5**2*S I GR < 2 - 5D0 > / 1 0 - 0 
GAM A=GAMA+DE X P < AL R05-ALR0 1 > *3 . 5**2*S I GR < 3 - 5D0 > / 1 0 . 0 
G AM A=G AM A+DE X P < A L ROG-ALRO 1 ) *4 . 5**2*S I GR ( 4 - 5D0 ) / 1 0 - 0 
G AM < I I ) =G AM A/ < P I + 1 97 . 0 > **2 

30 CONTINUE 
GAMAN=0.0 
PR0B=0.0 

DO 31 1 1=1, ISP 
GAMAG=GAMAG+GAM< I I >*RHO< I I ) 

PROB=PROB+RHOCII) 

31 CONTINUE 

GAM AG=G AM AG/ PROB 
GAMT=GA!iAG*3 - 2E21 
TAUG=1 .0/GAMT 
RETURN 
END 



